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A  reregulation  dam  Is  used  to  reduce  drastic  flow  fluctuations  from  main  reservoir 
releases,  usually  as  a  consequence  of  hydropower  operations.  In  a  pumped-storage  system, 
it  also  provides  storage  for  subsequent  pumpback  to  the  upstream  reservoir.  This  report 
describes  the  expected  impact  of  the  addition  of  a  reregulation  facility  on  the  tempera¬ 
ture  regimen  of  the  reservoir-reregulation  pool  system.  A  cost-effective  method  for  simu¬ 
lating  the  hydrodynamics  and  transport  in  a  reregulation  pool  is  presented.  The  method 
has  been  coupled  with  a  vertical  one-dimensional  reservoir  model  to  simulate  the  hydro¬ 
dynamics  and  transport  of  the  reservoir-reregulation  pool  system. 

The  modeling  procedure  for  simulating  the  hydrodynamics  and  transport  of  the  reregu¬ 
lation  pool  was  selected  from  three  candidate  methods.  The  three  methods  in  progressive 
order  of  complexity  are:  (1)  the  sump  method,  (2)  the  level  pool  routing  method,  and 
(3)  the  finite  difference  simulation  of  the  Saint-Venant  equations.  The  methods  have  been 
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incorporated  within  the  WESTEX  model,  and  the  predicted  temperature  profiles  in  the  main 
reservoir  and  the  release  temperatures  from  the  reregulation  pool  are  compared  to  results 
from  a  reservoir  study  previously  conducted  at  the  Waterways  Experiment  Station.  The 
longitudinal  resolution  of  the  level  pool  routing  and  the  Raint-Venant  simulation  methods 
allow  a  more  realistic  response  to  temperature  fluctuations  than  can  be  provided  by  the 
sump  method.  Further,  the  level  pool  routing  method  gives  results  comparable  to  the  more 
complex  Saint-Venant  method  in  a  more  cost-effective  manner. 


PREFACE 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  in  this  report 
(metric)  units  as  follows: 

_ Multiply _  _ _ 


acres 

acre-feet 

BTU  (International  Table) 
cubic  feet 

cubic  feet  per  second 
Fahrenheit  degrees 


4046.873 

1233.489 

1055.056 

0.02831685 

0.02831685 
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feet 

feet  per  mile 

*  2 

feet  per  second 

miles  (US  statute) 

pounds  (mass)  per  cubic 
foot 

square  feet 


0.3048 
0. 1893935 
9.290304 
1.609347 
16.01846 

0.09290304 


can  be  converted  to  SI 


_ To  Obtain _ 

square  metres 
cubic  metres 
joules 

cubic  metres 

cubic  metres  per  second 

Celsius  degrees  or 
Kelvins* 

metres 

metres  per  kilometre 
square  metres  per  second 
kilometres 
kilograms  per  cubic  metres 

square  miles 


*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F) 
readings,  use  the  following  formula:  C  =  (5/9) (F  -  32).  To  obtain 
Kelvin  (K)  readings,  use:  K  =  (5/9) (F  -  32)  +  273.15. 


HYDRODYNAMICS  AND  MODELING  OF  REREGULATION  POOLS 


PART  I:  INTRODUCTION 


Background 

1.  Gaining  the  full  potential  of  our  water  resources  can  often 
cause  conflicts  between  the  various  uses.  In  the  case  of  a  hydropower 
or  water-supply  reservoir,  the  demands  for  each  use  may  be  at  the 
expense  of  conditions  downstream.  For  example,  consider  Figure  1  in 
which  hydropower  operations  designed  for  peaking  power  production 
generally  induce  marked  pulsations  in  flow  during  the  day;  further, 


Figure  1.  Layout  of  reservoir-reregulatlon  pool  system  and 
the  smoothing  of  hydropower  operations 


» 

I 

I 

weekend  shutdowns  may  drastically  curtail  flow  in  the  downstream  area 
I  causing  wide  fluctuations  in  downstream  water  quality  (particularly  tern- 

i 

j  perature) .  This  wide  flow  variation  can  result  in  environmental  shock 

I  from  the  rapid  change  in  temperature  and  dissolved  oxygen.  Auxiliary 

i 

hydraulic  structures,  referred  to  as  reregulation  structures  (Figure  1), 
are  often  employed  just  downstream  of  hydropower  facilities  to  impound 
I  the  highly  variable  releases  from  these  projects.  Thus,  the  presence  of 

the  reregulation  pool  will  alleviate  these  downstream  problems  to  some 
degree . 

2.  In  addition  to  the  reregulation  of  flow,  some  reregulation 
pools  are  used  as  a  part  of  a  pumped-storage  system.  In  this  system, 
the  reregulation  pool  stores  water  from  generation  activities  during 
peak  demand  periods  so  that  some  portion  of  it  can  be  pumped  back  into 
the  upstream  reservoir  during  off-peak  hours.  In  this  manner,  the  water 
quality  of  each  reservoir  becomes  dependent  upon  the  other. 

3.  The  design  of  a  reregulation  pool  should  therefore  consider 
more  than  just  the  proper  storage  to  provide  a  uniform  discharge.  The 
effect  of  the  pool  upon  both  upstream  (for  pumped  storage)  and  down¬ 
stream  water  quality  should  be  analyzed.  Otherwise,  there  may  be  cause 
for  operational  revisions  which  could  impair  the  efficient  use  of  the 
reservoir  system. 

Objective 

4.  The  purpose  of  this  study  was  to  provide  insight  into  the 
potential  temperature  impact  of  a  reregulation  facility  and  to  develop 
an  effective  hydrodynamic  model  to  provide  input  in  the  simulation  of 
the  transport  of  temperature  and  conservative  constituents  in  a 
reservoir-reregulation  pool  system. 

Scope 


5.  The  potential  impact  of  a  reregulation  pool  on  water-quality 
constituents  will  be  demonstrated  through  the  use  of  previous 
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PART  II:  REVIEW  OF  I.T TFRATL'RE 


Reservoir  Dynamics 

6.  The  proper  simulation  of  the  reregulation  pool  and  Its  inter¬ 
action  with  the  surrounding  systems  requires  an  understanding  of  the 
major  physical  phenomena  affecting  the  reservoir-reregulation  pool 
system.  Such  understanding  will  allow  the  important  features  to  be 
simulated  at  the  resolution  necessary.  Thermal  stratification  is  a 
major  reservoir  feature;  it  has  a  profound  influence  upon  the  reser¬ 
voir's  circulation  and  the  resulting  water  quality  distribution.  Den¬ 
sity  stratification  is  a  result  of  the  effect  of  temperature  on  water 
density,  the  low  thermal  conductivity  of  water,  the  limited  penetration 
of  heat  through  the  air-water  interface,  and  the  seasonal  variation  of 
stream  temperatures  and  meteorological  conditions.  The  mechanics  of 
stratification  have  been  described  by  Hutchinson  (1957)  and  Kittrell 
(1965).  They  conclude  that  at  the  beginning  of  spring  the  reservoir  has 
a  virtually  uniform  vertical  density  profile.  During  spring  the  rate  at 
which  heat  is  distributed  downward  is  soon  exceeded  by  the  rate  of  in¬ 
coming  heat.  Most  of  the  energy  is  absorbed  near  the  surface.  Since 
warm  water  is  less  dense  it  floats  over  the  cool  water  and  a  definite 
temperature  gradient  develops,  thus  this  gradient  increases  during  the 
summer  months  of  high  energy  input.  The  diurnal  climatic  temperature 
fluctuations  and  wind  cause  mixing  of  the  surface  layers  forming  a  sur¬ 
face  region  of  relatively  uniform  temperature  called  the  epilimnion. 

The  deeper  regions  of  the  pool  remain  cool  due  to  the  stable  density 
stratification  which  has  developed.  This  lower  region  of  near  uniform 
density  is  called  the  hypolimnion.  The  transition  region  is  termed  the 
metalimnion.  Within  the  metalimnion,  the  elevation  which  has  the  great¬ 
est  vertical  temperature  gradient  is  called  the  thermocline.  In  the 
fall,  the  epilimnion  grows  as  a  result  of  surface  cooling  which  causes 
the  denser  surface  water  to  sink.  This  causes  mixing  which  eventually 
results  in  isothermal  conditions  again  by  the  beginning  of  spring. 


/.  Determination  of  whether  or  nor 
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for  particular  conditions  Is  somewhat  dlfiic’i’t. 
few  methods  that  give  some  general  i  r  d  1  ca  1 1  -n<-  . 
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fA  =  interfacial  friction  factor 

0  -  condenser  flow  rate 

c 

P  -  coefficient  of  thermal  expansion  for  water 

M  -  total  temperature  difference  fT  -  T  )  across  pond 
n  o  1 

T  =  discharge  temperature  (into  pond'* 

T  -  intake  temperature  (from  pond! 
g  -•  gravitational  acceleration  constant 
F  -  average  pond  depth 
W  =  average  pond  width 
D  =  vertical  entrance  dilutior 


v 

I.  =  total  pond  length 

Jirkn  and  Harleman  (19791  have  found  that  cooling  ponds  obey  a  one- 
dimensional  temperature  structure  (onlv  longitudinal  variation  exists! 
when  T./V  >  4  and  P  ;r  1.0  .  Ponds  for  which  1,/W  >  4  and  0.3  S  P  $ 
1.0  mav  he  analyzed  as  one  dimensional  (longitudinal)  if  allowance  is 
made  for  thermal  stratification  in  evaluating  the  surface-heat  flux. 
Ponds  in  which  P  c  n.3  were  found  to  he  vertically  well  stratified. 


*  For  convenience ,  svmbols  and  abbreviations  are  listed  In  the  Notation 
(Appendix  t ) )  . 


9.  Orlob  (1969)  suggested  the  use  of  the  densimetric  Froude 
number  to  assess  the  potential  for  vertical  stratification  in  a  lake. 
The  number  is  evaluated  as: 


where 

FD  ■  densimetric  Froude  number 
L  “  length  of  the  reservoir  in  ft 
Q  ■  flow-through  rate  in  cfs 
H  *  mean  reservoir  depth  in  ft 

3 

V  ■  volume  of  the  reservoir  in  ft 

2 

g  -  gravitational  acceleration  32.2  ft/sec  * 

e  -  average  normalized  density  gradient  taken  as  0.3  *  10  °/ft 
The  value  of  the  densimetric  Froude  number  that  marks  the  transition 
between  the  dominance  of  inertial  forces  and  the  dominance  of  gravita¬ 
tional  forces  Orlob  takes  to  be  1/ir  .  For  values  below  1/n  ,  a  strong 
vertical  stratification  pattern  will  exist. 

10.  Huber  and  Harleman  (1968)  listed  the  ratio  R  of  yearly 
reservoir  inflow  volume  to  actual  reservoir  volume  for  a  number  of 
impoundments  in  the  continental  United  States.  A  critical  value  R£  of 

about  5  to  10  has  been  suggested  so  that  R  <  R  could  be  considered  a 

c 

criterion  for  the  presence  of  reservoir  stratification  on  an  annual 
basis.  This  implies  that  a  reservoir  is  stratified  on  an  annual  basis 
only  if  the  detention  time  (volume  discharge)  is  greater  than  one  or  two 
months.  By  this  indicator  most  reregulation  pools  should  behave  in  a 
vertically  unstratified  manner  with  only  longitudinal  temperature  gradi¬ 
ents.  However,  the  main  reservoirs  are  generally  sufficiently  large  so 
that  an  annual  stratification  cycle  exists.  This  is  demonstrated  by 
Table  1  which  contains  a  compilation  of  US  Army  Corps  of  Engineers  data 

*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 


on  reservoir-reregulatlon  pool  systems.  Three  of  the  facilities  are 
used  In  conjunction  with  pumped-storage  operations.  The  reregulation 
pool  Is  generally  a  small  fraction  of  the  volume  of  the  main  reservoir 
for  each  project.  The  detention  times  for  the  reregulation  pools  in 
this  table  are  typically  less  than  20  days.  Therefore,  the  important 
gradients  are  in  the  vertical  direction  for  the  main  reservoir,  and  gen¬ 
erally  in  the  longitudinal  direction  in  the  reregulation  pool. 


Numerical  Description  of  the  Reservoir 

11.  In  order  to  study  the  impact  of  structural  or  operational 
modifications  on  the  reservoir  temperature  distribution  and  release,  nu¬ 
merical  models  which  reproduce  the  major  features  of  thermal  stratifica¬ 
tion  were  developed.  These  models  had  to  be  inexpensive  so  that  simula¬ 
tion  over  the  long  time  periods  required  to  reproduce  the  stratification 
cycle  of  the  reservoir  would  be  practical. 

12.  As  a  result,  the  US  Army  Engineer  Waterways  Experiment  Sta¬ 
tion  (WES)  developed  a  one-dimensional  thermal  model  (WESTEX,  based  pri¬ 
marily  upon  the  work  of  Clay  and  Fruh  (1970))  to  describe  the  stratifi¬ 
cation  process.  The  model  impoundment  is  discretized  into  horizontal 
layers  which  are  assumed  to  be  homogeneous.  A  particular  mass  of  water 
entering  a  stratified  reservoir  has  been  traced  with  dye  in  the  field 
(Elder  and  Wunderlich  1968,  Knight  1965)  and  found  to  eventually  reside 
at  the  vertical  location  with  approximately  the  same  density  as  the  in¬ 
flow.  Therefore,  the  WESTEX  model  places  the  inflow  in  the  layer  which 
is  closest  in  density.  The  velocity  profile,  due  to  outflow  from  the 
reservoir,  determines  the  contribution  from  each  layer  to  the  tempera¬ 
ture  of  the  release.  The  profile  is  determined  by  a  description  devel¬ 
oped  by  Bohan  and  Grace  (1973).  Edinger  and  Geyer  (1965)  and  Edinger, 
Brady,  and  Geyer  (1974)  developed  a  method  for  defining  the  heat  trans¬ 
fer  at  the  water  surface.  The  method  is  based  on  the  concept  of  an 
equilibrium  surface  temperature  towards  which  the  heat  exchange  process 
attempts  to  adjust.  This  has  been  incorporated  Into  WESTEX,  together 
with  an  exponential  absorption  of  short-wave  solar  radiation  with  depth 


due  to  light  penetration.  The  latter  is  based  upon  the  developments  of 
Dake  and  Harleman  (1966).  The  model  utilizes  an  extension  of  molecular 
diffusion  to  describe  the  nature  of  the  internal  heat  exchange  process. 
Internal  mixing  is  represented  by  a  mixing  scheme  based  on  an  integral 
energy  model  (Ford  1976). 

13.  Fontane  and  Bohan  (1974)  incorporated  a  second  reservoir 
capability  (for  simulating  the  reregulation  pool)  in  WESTEX  which  is 
essentially  a  zero  dimensional  thermal  budget  model.  This  reregulation 
pool  simulator  was  coupled  with  the  main  reservoir  through  the  main 
reservoir  releases  and  pumpback  flow  from  the  reregulation  model.  The 
model  system  continued  to  operate  inexpensively  since  the  modeling  pro¬ 
cedure  allowed  many  operations  (generation,  pumpback)  to  be  made  in  the 
reregulation  simulator  with  only  a  single  time  step  per  day  of  the  main 
reservoir  simulation.  This  reregulation  model  system  is  applicable  only 
where  the  assumption  of  homogeneity  (no  lateral,  longitudinal,  or  verti¬ 
cal  stratif ication)  of  the  reregulation  pool  is  reasonable.  This 
assumption  was  based  upon  physical  model  studies  and  field  data  that  in 
many  cases  indicated  that  the  reregulation  pool  was  not  vertically 
stratified  (Fontane  and  Bohan  1974,  Dortch  et  al.  1976,  Fontane  et  al. 
1977)  . 

14.  Not  only  does  the  residence  time  of  the  reregulation  pool 
affect  temperature,  but  the  release  flow  rate  from  the  pool  can  affect 
downstream  temperature.  Benedict  (1980)  modeled  the  Chattahoochee  River 
proposed  reregulation  pool  as  a  one-dimensional  longitudinal  model.  His 
model  indicated  that  reregulation  would  increase  the  summer  temperatures 
in  the  reach  of  the  pool  but  that  the  higher  base  flow  from  the  reregu¬ 
lation  pool  tended  to  offset  this  temperature  increase  downstream. 

15.  The  most  striking  result  of  a  reregulation  pool  is  the  impact 
upon  the  main  reservoir,  if  pumpback  is  used.  The  effect  of  pumpback 
from  the  reregulation  pool  upon  the  stratification  cycle  of  the  upper 
reservoir  has  been  recorded  in  numerous  studies.  Many  studies  have 
shown  that  pumpback  resulted  in  a  depressed  thermocline  (Simmons  1976; 
Oliver,  Hudson,  and  Clayton  1977;  Chen  and  Orlob  1972;  Potter,  Stevens, 
and  Meyer  1982).  The  results  of  studies  at  the  WES  (Holland  and  Dortch 


1982,  Fontane  et  al,  1977,  Dortch  1981)  also  indicated  that  a  more  homo¬ 
geneous  density  profile  results  in  the  main  reservoir,  following  the 
incorporation  of  pumped-storage  operations. 

16.  A  good  description  of  the  pumpback  entrainment  process  and  its 
resulting  effect  upon  stratification  is  given  by  Roberts  (1981).  Fig¬ 
ure  2  illustrates  the  phases  which  he  states  reduce  stratification.  A 
typical  summer  stratification  condition  in  the  reservoir  is  shown  in 
Figure  2a.  Following  pumpback  into  the  main  reservoir,  a  buoyant  jet 
(Figure  2b)  ascends  entraining  fluid  from  the  surrounding  reservoir. 

The  volume  flow  rate  of  the  jet  therefore  increases.  The  jet  then 
collapses  and  spreads  horizontally. 

17.  After  the  pumpback  process  ceases  (Figure  2c),  the  density 
current  formed  by  the  pumpback  jet  continues  to  spread.  The  resulting 
density  profile  will  be  less  stratified,  i.e.,  more  nearly  uniform  (Fig¬ 
ure  2d)  than  previously.  The  extent  of  the  impact  of  pumpback 


b.  Buoyant  jetting  inflow 


d.  Long  time  after  inflow 
ceases 


a.  Initial  conditions,  typical 
density  stratification 
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c.  Shortly  after  inflow  ceases 


Figure  2.  Hydrodynamic  features  of  pumped-storage  reservoirs 
subject  to  jetting  inflows  (from  Roberts  1981) 


entrainment  is  primarily  dependent  on  the  size  of  the  reservoir  (a  large 
reservoir  being  less  affected),  the  rate  of  the  pumpback  flow,  and  the 
strength  of  density  stratification  in  the  main  reservoir. 
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Numerical  Flow  Description  for  the  Reregulation  Pool 


18.  The  thermal  simulation  of  a  reregulation  pool  in  conjunction 
with  that  of  the  main  reservoir  needs  to  be  as  inexpensive  as  possible. 

It  would  appear  that  an  extension  of  the  WESTEX  code  for  coupling  a 
reregulation  pool  that  can  handle  longitudinal  gradients  requires  at 
least  a  one-dimensional  model.  The  precision  and  expense  of  the  simu¬ 
lation  method  hinges  upon  the  manner  in  which  the  hydrodynamics  (which 
can  involve  the  solution  of  two  nonlinear  partial  differential  equa¬ 
tions)  are  calculated.  If  the  density  driven  currents  are  considered 
small,  then  the  transport  of  heat  is  calculated  without  being  coupled 
with  the  hydrodynamics  and  involves  only  one  linear  partial  differential 
equation . 

19.  The  description  of  the  one-dimensional  equations  of  open  chan¬ 
nel  unsteady  flow  were  derived  by  Saint-Venant  (1871).  These  equations 
consist  of  the  conservation  of  mass  (continuity)  equation: 


3AV  3A 
3x  3t 
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and  the  conservation  of  momentum  equation  (equation  of  motion): 


3t 


(A) 


where 

A  ■  cross-sectional  area 
V  *  velocity 

x  *  distance  along  the  longitudinal  axis  of  the  waterway 


g  =  accelerlation  of  gravity 
h  =  water-surface  elevation  above  datum 


=  friction  slope 

20.  The  solution  of  the  set  of  partial  differential  equations 
requires  a  numerical  procedure  which  is  somewhat  laborious.  An  inter¬ 
mediate  approximation  uses  what  is  termed  hydrologic  modeling.  These 
models  are  based  on  the  conservation  of  mass  equation  in  the  following 
form  (Chow  1959): 
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in  which  AS  is  the  change  of  storage  within  the  reach  during  a  time 
increment;  the  storage  (S)  is  assumed  to  be  related  to  inflow  (I)  and/or 
outflow  (0) . 


S  =  K[XI  +  (1  -  X)0] 


(6) 


is  a 


K  is  a  storage  constant  with  dimensions  of  time  (secs),  and  X 
weighting  coefficient,  0  S  X  S  I  . 

21.  Storage  routing  models  attributed  to  Puls  (1928)  and  Goodrich 
(1931)  were  developed  assuming  X  to  be  zero.  Writing  Equation  3  in 
centered  time  and  space  finite  difference  form  yields  the  following 
reservoir  routing  model: 


I(t)  +  I(t  +  At) 
2 


0(t)  +  0(t  +  At) 
2 


S(t  +  At)  -  S (t) 
At 


(7) 


This  may  be  used  to  calculate  flow  distribution  throughout  the  pool 
assuming  a  level  pool  (Henderson  1966) . 

22.  Solving  the  actual  Saint-Venant  equations  has  been  addressed 
by  several  numerical  approaches,  one  of  the  earliest  being  the  charac¬ 
teristic  method.  The  basic  method  reduces  the  set  of  two  partial  dif¬ 
ferential  equations  to  a  set  of  four  ordinary  differential  equations  by 
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choosing  the  x  or  the  t  increment  so  that  the  Ax/At  ratio  lies 
along  a  characteristic  of  the  original  set  of  partial  differential  equa¬ 
tions.  Many  of  the  early  characteristic  schemes  were  explicit  (only  one 
unknown  per  equation),  such  as  the  schemes  of  Liggett  and  Woolhiser 
(1967),  Streeter  and  Wylie  (1967),  Lai  (1967),  and  Ellis  (1970). 

Implicit  (more  than  one  unknown  per  equation)  characteristic  models  were 
reported  by  Amein  (1966)  and  Wylie  (1970). 

23.  Many  finite  difference  numerical  models  of  the  Saint-Venant 
equations  have  been  reported.  In  finite  differences  the  partial  deriva¬ 
tives  are  replaced  by  difference  expressions  which  are  approximations. 
Early  work  was  registered  by  Stoker  (1953)  with  later  models  developed 
by  Garrison,  Granju,  and  Price  (1969);  Liggett  and  Woolhiser  (1967); 
Martin  and  De  Fazio  (1969);  Strelkoff  (3969);  Dronkers  (1969);  and 
Johnson  (1974).  There  are  of  course  many  variations  in  the  manner  in 
which  the  differences  are  formulated.  However,  all  the  methods  used  in 
the  references  above  suffer  stability  restrictions  which  result  in  small 
increments  of  t  .  The  limitation  is  roughly  that  a  surface  wave  cannot 
move  more  than  one  computational  cell  per  time  step.  As  an  example, 
consider  an  x  increment  of  1/2  mile  in  a  channel  30  ft  deep.  This 
criterion  would  limit  the  time  step  to  about  90  sec.  A  typical  reser¬ 
voir  simulation  period  is  1  year,  and  3.5  x  10^  time  steps  would  be 
needed.  Therefore,  explicit  schemes  would  not  be  practical  for  this 
study. 

24.  Implicit  finite  difference  schemes  solve  a  system  of  equa¬ 
tions  containing  more  than  one  unknown  per  equation  for  values  at  each 

X  grid  location  per  time  step.  Implicit  models  were  first  suggested  by 
Isaacson,  Stoker,  and  Troesch  (1956).  Many  authors  have  since  devised 
additional  implicit  schemes,  e.g.,  Preissmann  (1961),  Baltzer  and  Lai 
(1968),  Dronkers  (1969),  Amein  and  Fang  (1970),  Fread  (1973).  These 
methods  can  be  shown  to  be  unconditionally  stable  for  the  linear  advec- 
tion  equation. 
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which  is  an  analogy  for  the  Saint-Venant  equations. 

25.  The  nonlinear  partial  differential  equations  can  be  approxi¬ 
mated  as  a  set  of  finite  difference  equations  in  an  implicit  fashion 
with  the  following  approximations  suggested  by  Fread  (1973): 
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K  = 


0(k"+1  +  K“+J)  0  -  9>K  +  k;+1) 


(11) 


where 

K  -  any  function  or  variable 
Ax  and  At  *  space  and  time  increments 

i  and  „  *  the  mesh  x  location  and  time 
0  ■  time-weighting  factor 

The  temporal  difference  is  time-space  centered  and  the  spatial  deriva¬ 
tives  and  variable  representation  are  space-centered  between  i  and 
i+l  .  A  weighting  factor  of  0  “  1  yields  the  fully  implicit  scheme 

used  by  Baltzer  and  Lai  (1968).  A  weighting  factor  of  0  =  1/2 

produces  the  time-centered  "box"  scheme  used  by  Amein  (1966),  Amein  and 
Fang  (1970),  and  Contractor  and  Wiggert  (1972). 

26.  Fread  (1973)  has  shown  that  although  the  scheme  is  uncondi¬ 
tionally  stable  for  G  2  1/2  for  the  linear  advection  equation,  non¬ 
linear  stability  is  more  limited.  Nonlinear  stability  was  found  to 

require  0  >  0.55  .  Fread  found  that  the  closer  0  is  to  0.5  ,  the 
lower  the  numerical  distortion.  Fread  (1971)  also  developed  an  effi¬ 
cient  solution  algorithm  for  solution  of  the  linearized  system  of  equa¬ 
tions  taking  advantage  of  the  banded  coefficient  matrix. 


PART  III:  METHODS  OF  ANALYSIS 


Analytic  Analysis 

27.  In  order  to  evaluate  the  limitations  and  strengths  of  the 
candidate  simulation  methods  and  demonstrate  the  potential  effects  of 
reregulation  upon  the  transport  of  a  constituent,  a  progression  of 
simplified  models  will  be  considered.  In  these  examples,  temperature 
will  be  the  constituent  under  consideration;  application  to  a  conserva¬ 
tive  constituent  could  easily  be  made  in  the  same  fashion. 

28.  First,  consider  a  single,  well-mixed  reservoir  in  which  the 
inflow  matches  the  outflow. 


""COHO 


Qj  =  inflcw  and  outflow  rate 

Tj  =  instantaneous  temperature  of  the  reservoir  and  also  the 
outflow  temperature 

Tj  =  inflow  temperature 

Vj  «  volume  of  the  reservoir 

Sj  =  source  or  sink 

The  rate  of  heat  transfer  into  and  out  of  the  system  should  equal  the 
rate  of  change  of  heat  within  the  volume,  dV^T^/dt  .  Since  the  volume 
is  constant,  this  becomes  Vj(dT^/dt)  .  So,  Vj(dT^/dt)  =  heat  flow 
rate  in  -  heat  flow  rate  out  +  sources  -  sinks 
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A  second  reservoir  may  be  considered  which  is  fed  by  the  outflow  from 
reservoir  1.  This  second  reservoir  would  obviously  be  represented  as 


dT2  ^1  Q1T2  +  ^2 

dt  =  V2  T1  V2  ~  V2 


where  the  subscript  2  represents  the  second  reservoir.  If  the  lower 
reservoir  pumps  back  into  the  upper  reservoir,  this  would  be  a  heat  flow 
into  the  upper  and  lower  reservoirs.  In  order  to  keep  constant  reser¬ 
voir  volumes,  the  discharges  must  be  redefined.  Let  represent  the 

flow  from  reservoir  1,  Q ?  represent  the  pumpback  rate,  and  in  order  to 
balance  the  flows,  let  =  (Q^  -  Q7)  =  which  will  be  the  net  flow 

into  and  out  of  the  system.  The  differential  equation  for  the  first 
reservoir  would  then  be 


dT  Q  Q  Q  S 
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dt  v0  I  Vj  2  ]  -  V 


whereas,  the  second  reservoir  would  be 
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29.  Consider  a  two-reservoir  system  as  shown  in  Figure  3.  Reser¬ 
voir  i  serves  as  the  major  reservoir  and  reservoir  2  as  the  reregulation 
pool.  Assuming  both  reservoirs  to  be  homogeneous  (well-mixed),  the 
system  of  differential  equations  which  describes  the  heat  budget  of  each 
reservoir  is: 
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A  =  surface  area  of  reservoir  n 
n 

V  =  volume  of  reservoir  n 
n 

=  temperature  of  reservoir  n 
Qj  =  flow  rate  in  generation  (from  reservoir  1  to  2) 

=  pumpback  rate  (flow  rate  from  2  to  1) 

=  net  flow  rate  through  the  entire  system,  (equal  to  the 
inflow  to  reservoir  1,  and  also  the  release  from 
reservoir  2) 

Y  =  specific  weight  of  water 

c  =  specific  heat  of  water 
P 

n  =  1  or  2,  denoting  the  main  reservoir  and  reregulation  pool, 
respectively 

20.  Surface-heat  exchange  is  determined  by  an  approach  developed 
by  Edinger  and  Geyer  (1965)  and  Edinger,  Brady,  and  Geyer  (1974).  The 
thermal  equation  quantifying  the  net  surface-heat  exchange  (after  some 


linearization)  is: 


Hr  =  K(Tr  -  T) 


where 


1  =  net  rate  of  surface-heat  transfer 

s 

K  =  coefficient  of  surface-heat  exchange 
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Tj,  =  equilibrium  temperature 
T  =  water-surface  temperature 

Equilibrium  temperature  is  defined  as  that  temperature  at  which  the  net 
rate  of  heat  exchange  between  the  water  surface  and  the  atmosphere  is 
zero . 

HEAT  EXCHANGE 
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Figure  3.  Two-reservoir  system  with  putnpback 

31.  Assuming  that  and  Aj  ,  A^  are  constants,  then 

Qj  -  0^  +  Q3  .  The  system  can  then  be  considered  steady-state  for  flov 
for  certain  time  frames  such  as  on  a  daily  basis.  The  general  solution 
of  Equations  17  and  18  is: 
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=  initia]  temperature  of  reservoir  1 
=  initial  temperature  of  reservoir  2 

Homogeneous  reservoir 

system,  conservative  constituents 

32.  First,  consider  the  condition  in  which  there  is  no  surface- 
heat  exchange.  The  temperature  therefore  behaves  as  a  conservative 
constituent.  How  does  the  reregulation  pool  affect  the  transport  of  a 
constituent?  A  simple  example  will  be  used  to  illustrate  this  influ¬ 
ence.  The  example  system  contains  a  reregulation  pool  with  an  initial 
temperature,  T^  ,  equal  to  the  temperature  of  the  inflow  to  the  sys¬ 
tem,  Tj  .  The  major  reservoir  (reservoir  1)  has  an  initial  tempera¬ 
ture,  Tj  ,  1  deg  higher  than  the  initial  temperature  of  the  reregula¬ 
tion  pool.  If  the  initial  reregulation  pool  temperature  is  considered  a 
reference  temperature  (background  temperature),  then: 


T  =0 


These  temperatures  are  only  relative  temperatures  relating  the  major 
reservoir  and  the  reregulation  pool.  The  ratio  of  net  flow  Q-  to  the 


volume  of  the  major  reservoir  Is  Q^/^j  *  0,1/unit  time;  and  the 

volume  of  the  reregulation  pool  to  the  major  reservoir  V ^/V j  =0.1  . 
Figure  4  illustrates  the  transport  of  a  conservative  constituent 
(temperature  in  this  example)  through  the  two  reservoir  systems. 
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Figure  4.  Typical  temperatvire  of  a  two- 
reservoir  system  for  Q^/V^  =  ^.1  , 

V  /V  =  o.l  ,  with  varying  Q ^/V ]  ratios 
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33.  The  example  shows  the  temperature  time  history  for  the  two 
reservoirs  for  three  levels  of  pumpback;  Q2 /V ^  =  0  ,  Q£ /V ^  =0.1  , 
and  (^/V^  =  *  •  The  temperature  of  the  major  reservoir  is,  of  course, 
at  a  maximum  initially  and  is  diluted  continuously  thereafter.  The  re¬ 
regulation  pool  temperature  would  rise  until  it  reaches  a  peak  when  it 
is  identical  to  the  major  reservoir's  temperature  at  that  time,  after 
which  it  also  declines.  As  is  illustrated,  the  peak  reregulation  pool 
temperature  occurs  more  quickly  as  the  amount  of  pumpback  is  increased. 
For  no  pumpback,  the  time  to  peak  (obtained  by  setting  the  time-rate 
derivative  of  Equation  21  to  zero)  is  given  by: 

In  (V  /V  ) 

Cp  =  D  Q1(V1  -  V2)/VLV2  (22 


As  pumpback  increases,  the  time  to  peak  is  reduced,  limited  by  t  =  0 


of  course.  In  the  example,  t  =2.56  for  the  no  pumpback  case;  1.54 

P 


for  C^/V^  =0.1;  and  0.41  for  (^/V ^  =  1  .  The  actual  peak  tempera¬ 


ture  in  the  reregulation  pool,  for  this  example,  increases  with  increas 
ing  pumpback.  In  this  example,  the  relative  peak  temperature  of  the 
reregulation  pool  was  0.76  for  no  pumpback;  0.82  for  (^/V^  —0.1  ;  an 


0.88  for  Q2/Vj 


1 


The  peak  temperature  approaches  (V^/(Vj  + 


Vr 2 ) ) T 1 1  as  pumpback  is  increased.  The  larger  the  pumpback  flow  rate, 
the  more  the  two  reservoirs  behave  as  if  they  were  a  single  reservoir 
which  has  a  volume  equal  to  the  sum  of  the  two  reservoirs  in  the  system 
34.  Figure  5  shows  the  peak  relative  temperature  and  the  time  to 
this  peak  in  the  reregulation  pool  as  a  function  of  the  ratio  of 


reregulation  pool  volume  to  major  reservoir  volume,  V ^/V 1  ,  for 


Q^/Vj  =0.1.  As  the  relative  volume  of  the  reregulation  pool  is  in¬ 
creased,  the  temperature  peak  is  attenuated.  The  time  to  peak  temper¬ 
ature  increases  sharply  as  the  reregulation  pool  is  enlarged.  It  is 
evident  from  Figure  5  that  pumpback  has  negligible  impact  on  the  tem¬ 
perature  in  the  system  if  the  reregulation  pool  is  small.  Obviously, 
this  applies  only  in  a  well-mixed  system  (no  stratification). 
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Figure  5.  Conservative  constituent,  Q~/Vj  =  0.1 

Homogeneous  reservoir  system, 
nonconservative  constituents 

35.  Figure  6  demonstrates  the  effect  of  adding  surface-heat 
exchange.  For  this  example  the  following  values  of  parameters  were 
used : 


Volume 

Surface 

Area 

Inflow  and 
Release, 

Pumpback, 


Main  Reservoir 
300  KAF* 

3  KA** 

0.5  KAF/Dt 

10.0  KAF/D 


Reregulation  Pool 
10  KAF 
1  KA 


*  KAF  =  kiloacre-f eet . 

**  KA  =  kiloacre. 
t  KAF/D  =  kiloacre-feet/day . 

The  surface-heat  exchange  coefficient  used  in  this  example  was  K  =  140 
? 

BTU/ft  -day-°F.  The  inflow  temperature  and  the  initial  temperaturi 
the  reservoirs  were  arbitrarily  set  at  1°  F  below  the  equilibrium 
temperature.  This  allowed  the  computation  of  H  as  discussed  ,-n 
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Figure  6.  Pumpback  effects  with  heat  exchange 

paragraphs  27  through  31.  In  the  absence  of  pumpback,  the  major 
reservoir  is  unaffected  by  the  reregulation  pool.  While  the  detention 
time  in  the  main  reservoir  is  quite  large  (600  days),  its  average  depth 
is  also  large  (100  ft),  thus  the  main  reservoir  heats  slowly.  The 
ultimate  steady-state  temperature  for  the  main  reservoir  (T^g  )  under 
these  conditions  (from  Equation  20  with  t  =  “>  )  would  be  given  by 
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or  0.930.  The  rate  constant  associated  with  the  upper  reservoir  is 

M  *  [Q  +  K/(vc  )A  ] / V  =  0.0241  .  With  no  pumpback,  the  rate  indica- 
i  1  Dll 


tors  are  Mj  =  -  [Qj  +  (K/ycp) ] /V ^  and  M2  *  -  [Qj  +  (K/ycp) A^ ] /V^  . 


The  magnitudes  of  M^  and  M2  are  the  rate  constants  associated  with 
the  reregulation  pool  and  the  main  reservoir,  respectively.  The  larger 
the  rate  constant,  the  more  rapidly  the  temperature  changes. 

36.  The  rate  constant  associated  with  the  reregulation  pool  is 
substantially  larger  than  that  of  the  main  reservoir  (0.274  compared 
with  0.0241),  so  the  reregulation  pool  is  more  rapidly  influenced  by  the 
equilibrium  temperature  than  the  main  reservoir.  Though  the  reregula¬ 
tion  pool  cannot  influence  the  main  reservoir  without  pumpback,  the 


reregulation  pool  is  fed  by  the  upper  reservoir  so  that  temperature 
changes  which  are  produced  upstream  are  built  upon  by  the  reregulation 
pool.  The  ultimate  steady-state  temperature  of  the  reregulation  pool 
(T2sg)  would  be 


For  this  example,  the  ?2Ss  =  ®*^'7  *  The  rate  of  change  of  temperature 
in  the  reregulation  pool  is  quite  rapid  compared  to  the  upstream  reser¬ 
voir.  Although  the  detention  time  in  the  reregulation  pool  is  only 
20  days  compared  to  600  in  the  upstream  reservoir,  the  average  depth  of 
the  reregulation  pool  is  one  tenth  that  of  the  main  reservoir;  this,  in 

turn,  allows  the  reregulation  pool  to  have  quicker  temperature  fluctua¬ 

tions  than  the  main  reservoir. 

37.  Figure  6  also  shows  the  effects  of  pumpback  at  a  rate  of 

10  KAF/D  ,  where  10  KAF  is  the  volume  of  the  reregulation  pool.  This 

rate  of  pumpback  resulted  in  a  substantially  different  reregulation  pool 
temperature  history  than  that  computed  for  no  pumpback.  The  impact  of 
this  pumpback  rate  on  the  main  reservoir  was,  however,  considerably  less 
than  that  on  the  reregulation  pool  due  to  the  large  size  of  the  main 
reservoir.  Again,  the  pumpback  caused  mixing  of  the  two  reservoirs  so 
that  they  began  to  behave  as  a  single  reservoir.  The  temperature  his¬ 
tory  for  a  pumpback  rate  =  00  is  also  shown  in  this  figure.  Under 
this  condition,  T^  =  T^  .  Since  the  reregulation  pool  is  much  smaller 
than  the  main  reservoir,  the  effect  of  pumpback  is  more  pronounced  upon 
it  than  upon  the  larger  reservoir.  An  important  feature  of  this  pump¬ 
back  is  that  during  periods  of  warming,  the  reregulation  pool  will 
generally  warm  more  quickly  than  the  main  reservoir,  and  the  pumpback 
will  tend  to  heat  the  main  reservoir. 

38.  The  presence  of  the  reregulation  pool  obviously  promotes  con¬ 
ditions  different  from  natural  conditions  within  this  same  reach.  The 


reregulation  dam  stores  pulses  of  flow  and  releases  a  relatively  uniform 
flow.  This  storage  produces  a  different  response  in  the  reach  to  the 
equilibrium  temperature  fluctuations.  Again  considering  daily  average 
events,  analysis  of  these  effects  would  be  quasi-steady  state.  The  sys¬ 
tem  may  be  described  by: 
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where 

T  ■  temperature  of  the  reach 
Q  ”  flow  through  the  reach 
V  “  volume  of  the  reach 
A  =  surface  area  of  the  reach 
a  sin  0t  “  equilibrium  temperature  as  a  function  of  time 
0  =  cycle  frequency  ■  2ir/T 

c 

T^  =  cycle  period 

The  solution  for  the  case  in  which  the  initial  temperature  of  the  pool 
is  assumed  to  be  zero  (compared  to  a  given  reference  point)  is: 


Set 
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where 
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39.  The  characteristics  of  a  reach  will  determine  the  system's 

temperature  response  to  the  particular  heating  cycle.  If  the  frequency 

of  the  heating  cycle  is  very  low  (i.e.f  1  year),  the  amplitude  ratio  of 

the  temperature  fluctuation  in  the  pool  to  that  of  the  environment  would 

then  be  given  by  (KA/yc  )/[Q  +  (KA/yc  )]  .  So,  if  Q  is  small 
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compared  to  the  heat  exchange  component  (KA/yc  ),  the  detention  time  in 

P 

the  reach  is  long  enough  for  the  range  of  fluctuations  of  the  temper¬ 
ature  in  the  reach  to  almost  match  that  of  the  environment.  Conversely, 
a  larger  discharge  through  the  reach  would  shrink  the  temperature  range. 
Thus,  the  increased  surface  area  associated  with  the  reregulation  pool 
tends  to  cause  greater  temperature  fluctuations  compared  to  base  (no 
reregulation  pool)  conditions  and  in  response  to  low  frequency  heating 
cycles.  For  high  frequency  loadings,  the  temperature  response  would  be 
quite  small  except  for  very  shallow  pools. 

40.  As  an  example,  consider  the  previous  instance  where  the  pool 
volume  was  10  KAF,  the  surface  area  was  1  KA,  the  discharge  was 

? 

0.5  KAF/D,  and  the  surface-heat  exchange  coefficient  was  140  BTU/ft^- 
day-°F.  Table  2  shows  the  ratio  of  the  amplitude  of  temperature  fluctu¬ 
ation  cycle  experienced  in  the  reregulation  pool  reach  to  that  of  the 
equilibrium  temperature  for  conditions  with  and  without  (natural)  a 
reregulation  pool.  The  atmospheric  heating  cycles  which  were  of  primary 
interest  were  1  day  and  1  year.  These  would  probably  describe  the 
strongest  cyclic  temperature  variations  for  the  reach.  This  table  fur¬ 
nishes  information  for  these  two  cycles  and  for  a  frequency  of  zero 
(infinite  heating  period).  The  critical  factor  for  the  high  frequency 
diurnal  heat  loading  (the  1-day  loading  period)  is  the  mean  depth  of  the 
pool  (in  this  example  10  ft).  Of  course,  a  small  mean  depth  implies  a 
large  surface  area  for  a  given  pool  volume.  Even  for  such  a  shallow 
depth,  the  response  of  the  pool  was  less  than  4  percent  of  the  equili¬ 
brium  temperature  amplitude.  The  annual  cycle  produced  a  response  in 
the  pool  of  81  percent  of  the  amplitude  of  the  equilibrium  temperature. 
So  for  all  but  the  shallowest  conditions,  the  daily  fluctuations  in  pool 
temperatures  would  be  quite  small  compared  to  that  of  the  annual  cycle 
for  the  reregulation  pool  reach. 

41.  As  shown  in  Figure  7,  a  schematized  natural  flow  condition 
would  be  within  the  channel  banks,  whereas  the  reregulation  pool  might 
often  flood  areas  above  this  elevation  and  spill  into  the  natural  flood 
plain.  For  example,  let  the  natural  surface  area  be  0.25  KA  and  the 
volume  be  5  KAF.  This  means  the  average  depth  is  20  ft.  Compared  to  a 
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Figure  7.  Schematized  cross  sec¬ 
tion  of  a  reregulation  pool 
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condition  with  a  flooded  overbank,  this  greater  depth  leads  to  even  less 
sensitivity  to  higher  frequency  heating  cycles  than  computed  for  the 
reregulation  pool.  Column  2  of  Table  2  shows  that  for  the  diurnal  cycle 
the  example  natural  pool  has  a  temperature  range  of  less  than  2  percent 
of  that  of  the  environment.  The  combination  of  increased  mean  depth  and 
a  reduction  in  detention  time  (resulting  from  all  flow  staying  in-bank 
for  the  natural  reach)  caused  the  temperature  amplitude  in  the  natural 
reach  to  be  only  53  percent  of  the  long-period  equilibrium  temperature 
cycle  amplitude.  The  reregulation  pool  produced  a  temperature  range 
ratio  of  over  80  percent  for  this  condition.  Therefore,  the  reach  in 
which  the  example  reregulation  pool  is  located  would  tend  to  have 
larger-than-natural  annual  and  diurnal  temperature  variations  due  to  its 
increased  detention  time  and  surface  area. 

42.  ViThile  the  reach  which  contains  the  reregulation  pool  may  have 
a  greater  temperature  fluctuation  than  occurs  without  the  reregulation 
pool,  this  is  not  necessarily  the  case  below  the  reregulation  dam.  The 
primary  function  of  the  dam  is  to  store  widely  varying  generation  flow 
and  produce  a  relatively  uniform  outflow.  The  extreme  low  flow  condi¬ 
tions  which  may  occur  without  the  reregulation  pool  during  nongeneration 
periods  could  allow  quicker  adjustment  of  the  water  temperatures  to  the 
equilibrium  temperature  environment  than  would  the  uniform  release  pro¬ 
duced  by  the  reregulation  pool.  The  water  would  become  more  acclimated 
to  the  equilibrium  temperature  conditions  during  nongeneration  periods 
due  to  the  shallower  mean  depth  and  decreased  velocities  than  those  for 
analogous  periods  of  higher  uniform  releases  from  the  reregulation  pool. 
Conversely,  during  the  periods  of  generation  activity  with  no  reregu¬ 
lation  pool,  the  flow  would  be  less  acclimated  to  the  temperature  of  the 
environment  than  would  be  experienced  with  the  lower,  albeit  uniform 
reregulation  pool  release.  Thus,  in  summer  or  warm  season  months  the 


temperature  range  experienced  in  the  reach  below  a  hydropower  facility 
with  no  reregulation  pool  could  be  quite  large  due  to  the  wide  variation 
in  flow  between  generation  and  nongeneration  periods.  This  variation 
would  be  lessened  somewhat  by  the  uniform  nature  of  reregulation  pool 
outflow. 

43.  These  conclusions  are  similar  to  the  results  described  by 
Benedict  (1980)  in  a  modeling  effort  concerning  the  proposed  Chattahoo¬ 
chee  reregulation  pool  below  Buford  Dam.  Benedict  was  concerned  with  a 
summer  heating  scenario  in  which  the  proposed  reregulation  pool,  as 
modeled,  would  increase  temperatures  in  the  reach  relative  to  existing 
conditions.  However,  flow  variations  in  the  reach  would  be  dampened  by 
the  reregulation  dam  so  the  temperature  fluctuations  below  the  dam  would 
be  reduced  relative  to  existing  conditions.  His  results  showed  that  the 
higher  base  flows  provided  by  the  reregulation  dam  would  tend  to  offset 
some  of  the  warming  caused  by  the  reregulation  pool.  Zimmerman  and 
Dortch  (1986)  modeled  the  Chattahoochee  River  in  an  unsteady  manner. 

They  found  that  modeling  of  the  dynamics  was  essential  to  register  down¬ 
stream  temperature  features  in  addition  to  those  due  to  higher  base 
flows. 

44.  The  discussion  thus  far  has  concerned  homogeneous  reservoirs. 
However,  potentially  the  most  significant  effect  of  a  reregulation  pool 
in  a  pumped-storage  hydropower  system  will  be  its  effect  upon  the 
stratification  cycle  of  the  main  reservoir.  As  was  previously  demon¬ 
strated,  a  large  reservoir  with  a  relatively  small  flow-through  has  a 
large  detention  time  and  would  have  a  wide  temperature  variation  for 
long-period  heating  cycles.  This  analysis  indicates  the  entire  reser¬ 
voir  should  get  quite  warm  in  the  summer  and  cool  in  the  winter.  How¬ 
ever,  this  is  not  generally  the  case.  In  the  spring  or  early  summer, 
heat  is  transferred  into  the  surface  layers  more  rapidly  than  it  can  be 
diffused  or  mixed  into  the  subsequent  lower  layers.  Since  warm  water  is 
less  dense  than  cold  water,  it  tends  to  rise  toward  the  surface.  This 
action  forms  a  density  gradient  which  effectively  inhibits  any  further 
mixing.  This  process  results  in  a  density-stratified  reservoir. 
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45.  As  the  cooling  season  approaches,  the  upper  layers  are  cooled 
and  thus  sink.  This  action  reduces  density  gradients  so  that  eventually 
circulation  is  often  produced  over  the  entire  depth  during  the  winter; 
the  reservoir  subsequently  becomes  relatively  homogeneous.  The  pumpback 
into  the  main  reservoir  can  reduce  the  stratification  in  the  upper 
reservoir  by  both  direct  and  indirect  means.  The  direct  manner  is  due 
to  the  mixing  and  entrainment  by  the  pumpback  plume  in  the  main  reser¬ 
voir.  The  indirect  method  is  due  to  the  return  of  the  warmer  water  from 
the  reregulation  pool  causing  a  further  reduction  of  the  density  gradi¬ 
ent.  This  return  flow  from  the  reregulation  pool  can  drastically  modify 
the  stratification  cycle  of  the  main  reservoir.  Thus,  the  effect  of 
stratification  upon  the  reregulation  pool  dynamics  must  be  considered. 
The  process  is  sufficiently  complex  so  that  its  analysis  usually 
requires  modeling. 


Summary  of  reregulation  pool 
effects  from  analytic  analysis 

46.  The  storage  associated  with  a  reregulation  pool  will  result 
in  longer  detention  rimes  than  would  occur  without  a  reregulation  pool. 
Therefore,  the  annual-cycle  temperature  amplitude  of  waters  in  the  reach 
could  be  increased  compared  to  nonrereregulated  conditions  and  result  in 
warmer  summer  water  temperatures  in  the  reach  of  the  pool. 

47.  Downstream  of  the  reregulation  pool,  the  more  uniform  flows 
produced  by  this  pool  may  moderate  short  period  (diurnal)  temperature 
fluctuations  compared  to  nonreregulated  conditions.  The  longer-period 
temperature  cycles  may  also  be  moderated  some  distance  below  the  reregu¬ 
lation  dam. 

48.  The  response  to  daily  temperature  cycles  in  the  reregulation 
pool  reach  will  generally  be  slight.  The  response  is  primarily  depen¬ 
dent  on  the  mean  depth;  therefore,  the  diurnal  fluctuation  in  tempera¬ 
ture  within  the  reregulation  pool  may  be  more  or  less  than  the  natural 
fluctuation,  depending  upon  the  bathymetry  of  the  pool. 

49.  In  a  pumped-storage  system,  large  pumpback  rates  cause  the 
upper  reservoir  and  reregulation  pool  to  become  similar  in  quality  if 
the  upper  reservoir  and  the  reregulation  pool  are  homogeneous.  During 
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warming  trends  the  reregulation  pool,  in  conjunction  with  pumpback 
operations,  will  tend  to  warm  the  main  reservoir  more  quickly.  The  con¬ 
centration  of  a  conservative  constituent  in  the  main  reservoir  will 
attain  the  same  value  in  the  reregulation  pool  more  quickly  with 
increasing  use  of  pumpback  generation.  The  peak  concentration  in  the 
reregulation  pool  will  again  tend  to  be  higher  with  greater  pumpback  and 
general ion . 

50.  Pumpback  into  the  main  reservoir  from  the  reregulation  pool 
can  result  in  considerably  more  uniform  density  profiles  there.  The 
degree  of  influence  of  the  entrainment  and  mixing  depend  on  the  relative 
size  of  the  upper  reservoir,  the  magnitude  of  pumpback,  and  the  strength 
of  density  stratification  in  the  main  reservoir. 

Numerical  Modeling  Analysis 

51.  The  interchange  of  water  between  the  upper  reservoir  and  the 
reregulation  pool  can  have  a  variety  of  effects  upon  the  environmental 
regimen  of  either  pool.  As  Clugston  and  Prince  (1980)  stated:  "The 
severity  of  the  alterations  depends  on  the  surface  area,  depth,  and 
morphometry  of  both  reservoirs;  the  volumes  of  water  exchanged;  water 
retention  time;  depth  of  intake  and  discharge  structures;  and  season  and 
time  of  day  during  which  a  plant  operates."  Consequently,  a  modeling 
system  that  has  sufficient  flexibility  and  is  relatively  inexpensive  for 
simulation  of  long  periods  of  interest  is  necessary  to  obtain  reliable 
results  for  actual  investigations.  This  particular  study  builds  upon 
previous  work  conducted  at  WES  in  which  the  thermal  model  WESTEX  was 
developed  and  extends  the  modeling  capability  associated  with  simulation 
of  the  reregulation  pool  in  tandem  with  the  main  reservoir. 

52.  The  WESTEX  model  is  a  one-dimensional  numerical  model  which 
provides  the  vertical  stratification  structure  needed  to  describe  hydro- 
dynamic  behavior  in  the  main  reservoir.  This  particular  modeling  effort 
will  concentrate  upon  the  incorporation  into  the  WESTEX  model  of  a 
numerical  one-dimensional  (longitudinal)  description  of  the  reregulation 
pool  to  provide  information  on  longitudinal  gradients.  A  basic 
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assumption  is  that  the  pool  is  well-mixed  vertically  and  laterally. 

This  assumption  seems  reasonable,  given  the  short  detention  times  asso¬ 
ciated  with  existing  reregulation  pools.  By  modeling  each  reservoir  as 
a  one-dimensional  system  (longitudinal  resolution  in  the  reregulation 
pool  and  vertical  in  the  main  reservoir),  the  physically  significant 
characteristics  of  each  reservoir  will  be  maintained  while  keeping  com¬ 
puter  costs  low. 

Discussion  of  WESTEX  model 

53.  The  basic  equation  which  is  solved  by  WESTEX  is  the  one¬ 
dimensional  advection-dif fusion  equation.  It  relates  the  increase  in 
heat  in  an  element  to  the  transport  of  heat  by  advection  or  diffusion 
across  the  boundaries  of  the  element  along  with  the  sources  and  sinks. 
This  equation  is  derived  as  follows: 

Advection : 


The  increase  in  heat  due  to  vertical  advection  is 

/  30  \  -3Q  T 

Atly  -  v  -  3r TAz; =  — AtAz 

where 

At  =  the  time  step  length 

Q  =  the  net  vertical  flow  into  or  out  of  a  layer 
v 

The  increase  in  heat  due  to  inflow  and  withdrawal  is 


where 

=  inflow  rate,  cfs/day 

=  temperature  of  the  inflow,  °F 

Q  =  withdrawal  rate,  cfs/day 
o 

T  =  temperature  of  the  withdrawal,  °F 
o 

Dispersion : 

The  heat  transfer  is  assumed  to  be  described  by  Fickian  diffusion. 
Transport  through  the  top  of  the  element  is  given  by 
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Whereas,  transport  through  the  bottom  of  the  element  is 
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Therefore,  the  heat  transfer  due  to  dispersion  can  be  written  as 

h  (“  i)4Z4t 

where 

2 

K  =  vertical  diffusion  coefficient,  ft  /day 

r  2 

A  =  horizontal  area  of  the  element,  ft 
External  source: 

The  increase  in  heat  due  to  an  external  source  is  given  by 
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where 
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p  =  density  of  water,  lb/ft 

c  =  specific  heat  of  water,  RTU/lb/°F 
P 


Vv 
V.'-J 

m 


34 


external  heat  source,  BTU/ft/day 


3Z 

54.  The  sum  of  these  terms  should  be  the  rate  of  increase  of 
temperature  multiplied  by  the  volume  of  the  element  control  volume 
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which  then  yields  the  following  equation 
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Now  dividing  by  AAZAt  results  in 
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This  expression  for  conservation  of  heat  is  coupled  with  the  conserva¬ 
tion  of  mass  equation  for  an  elemental  control  volume.  The  model  code 
was  developed  with  consideration  for  the  long  time  periods  required  to 
simulate  the  natural  stratification  cycle  of  a  reservoir.  For  this 
reason,  a  one-dimensional  approach  was  followed.  WESTEX  was  developed 
at  WES  but  was  based  upon  the  prior  work  of  Bohan  and  Grace  (1969, 
1973);  Clay  and  Fruh  (1970);  Edinger  and  Geyer  (1965);  Edinger,  Brady, 
and  Geyer  (1974);  and  Dake  and  Harleman  (1966).  The  model  includes 
computational  methods  for  simulating  heat  transfer  at  the  air-water 
interface,  advective  heat  transfer  of  inflow  and  outflow,  and  internal 
dispersion  of  thermal  energy.  Since  the  model  is  conceptually  one¬ 
dimensional,  the  impoundment  is  divided  into  discrete  horizontal  layers 
of  uniform  thickness.  Fundamental  assumptions  include  the  following: 

a.  Isotherms  are  parallel  to  the  water  surface. 

b.  The  water  in  each  discrete  layer  is  isotropic  and 
physically  homogeneous. 
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d.  External  advection  (inflow  and  outflow)  occurs  as  a 
uniform  distribution  within  each  layer. 

£.  Internal  dispersion  of  thermal  energy  is  accomplished  by 
a  lumped  diffusion  mechanism  that  combines  the  effects 
of  molecular  diffusion,  turbulent  diffusion,  and  thermal 
convection. 

£.  There  is  no  heat  exchange  through  the  mud-water  interface 
at  the  bottom. 

A  description  of  this  model  can  be  found  in  Holland  and  Dortch  (1982). 

55.  The  WESTEX  model  has  used  a  sump  model  to  simulate  a  re¬ 
regulation  pool.  This  procedure  for  reregulation  pool  (or  afterbay) 
simulation  is  the  approach  used  in  past  applications  of  WESTEX.  Its  use 
is  documented  in  several  previous  investigations  (Fontane  et  al.  1977, 
Dortch  et  al.  1976,  and  Holland  and  Dortch  1982).  The  reregulation  pool 
is  modeled  numerically  by  maintaining  heat  and  water  budgets.  Pumpback 
would  logically  retrieve  some  of  the  water  recently  entering  the  after¬ 
bay  during  a  generation  cycle.  This  method  utilizes  a  pumpback  coeffi¬ 
cient  which  is  the  percentage  of  the  pumpback  water  volume  that  is  from 
the  day's  generation.  The  rest  of  the  pumpback  amount  comes  from  the 
afterbay.  The  remainder  of  generation  volume  is  then  fully  mixed  with 
the  afterbay.  In  this  manner,  some  allowance  is  made  to  account  for 
longitudinal  gradients.  However,  a  single  sump  system  such  as  this  has 
only  a  one-day  memory  since  the  reregulation  pool  is  fully  mixed  at  the 
end  of  each  day.  The  system  requires  some  estimation  (through  a  coeffi¬ 
cient)  of  the  percentage  of  generation  flow  making  up  the  pumpback 
volume.  The  two  extremes  in  this  estimation  are:  (a)  plug  flow  (no 
longitudinal  mixing),  where  the  pumpback  volume  is  made  entirely  of 
generated  volume  and  (b)  fully  mixed,  where  the  pumphack  volume  is  made 
up  of  well-mixed  reregulation  pool  quality  after  generation 

56.  Once  during  each  day  of  simulation  the  volume  of  water  in  the 
reregulation  pool  is  adjusted  to  sequentially  account  for  (a)  the  main 
reservoir  generation  volume  which  is  not  pumped  back;  (b)  the  afterbay 
volume  which  is  pumped  back;  and  (c)  the  reregulation  pool  volume  which 
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is  released  downstream.  The  net  contributions  of  the  generation  volumes 
from  the  main  reservoir  for  each  operation  cycle  at  their  respective 
temperatures  are  added  to  the  post  pumpback  afterbay  volume  and  a  volume- 
weight  average  temperature  for  the  afterbay  is  computed.  The  daily 
surface-heat  exchange  is  then  applied  to  the  afterbay  surface  area  and  a 
new  afterbay  temperature  computed.  The  temperature  of  the  pumpback 
volume  is  computed  as  a  flow-weighted  average  of  the  mixed  temperatures 
of  the  generation  volumes  pumped  back  and  the  percentage  of  the  reregu¬ 
lation  pool  (at  its  temperature)  which  is  pumped  back. 

57.  An  example  of  this  procedure  may  illustrate  the  method  more 
clearly.  Consider  for  example: 

a.  A  reregulation  pool  with  a  volume  of  10  KAF  and  a 
temperature  of  60°  F  at  the  beginning  of  the  day. 

b.  A  generation  rate  of  10  KAF/D  for  one-half  day  which 
withdraws  water  of  40°  F  temperature  from  the  main 
reservoir. 

c.  A  pumpback  rate  of  8  KAF/D  for  one-half  day  with  a  pump¬ 
back  coefficient  of  0.6. 

d.  The  release  rate  from  the  reregulation  pool  is  1  KAF/D 
for  the  entire  day. 

The  volumes  added  or  removed  from  the  reregulation  pool  are  simply  the 
flow  rate  multiplied  by  the  duration;  for  this  example,  5  KAF  for 
generation,  4  KAF  for  pumpback,  and  1  KAF  for  the  release  downstream. 

The  temperature  of  the  pumpback  is  found  by  the  flow-weighted  average 
temperature  of  the  generation  flow  returned  to  the  main  reservoir  (de¬ 
termined  by  the  pumpback  coefficient),  and  the  remainder  of  the  pumpback 
flow  which  is  from  the  reregulation  pool.  In  this  example,  the  pumpback 
temperature  is  given  by: 

Pumpback  temperature  =  (0.6) (40°  F)  +  (0.4) (60°  F)  =  48°  F 

Of  course,  the  maximum  amount  of  generation  flow  which  can  be  returned 
is  limited  to  that  day's  generation  volume.  The  temperature  and  volume 
of  the  reregulation  pool  is  then  calculated  by  summation  of  the  flow 
volume  and  the  flow-weighted  temperature  of  each  component  of  the 
activities  of  the  day.  In  this  example, 


Reregulation  pool  volume  =(10+5-4-1)  KAF  =  10  KAF 


Reregulation  pool  volume  =  10  KAF  +  5  KAF  -  4  KAF  -  1  KAF  =  10  KAF 


Reregulation  pool  temperature 

_  [ ( 10KAF) (60°  F)  +  (5KAF) (40°  F)  -  (4KAF)(48°  F)  -  (1KAF)(60°  F) ] 

10  KAF 

=  54.8°  F 

At  this  point  the  surface-heat  exchange  for  the  day  is  applied. 

Methods  for  simulating 
reregulation  pools  hydrodynamics 

58.  An  alternative  to  the  sump  model  which  provides  longitudinal 
resolution  calculates  the  current  velocity  and  depths  by  level  pool 
routing  and  the  constituent  transport  is  then  calculated  by  a  one¬ 
dimensional  transport  model.  As  the  name  implies,  in  the  calculation  of 
hydrodynamic  parameters  this  modeling  concept  assumes  the  pool  water 
surface  to  be  level.  With  this  assumption,  only  the  lumped,  unsteady 
continuity  equation  has  to  be  solved  to  compute  discharge  at  each  cross 
section.  The  method  proceeds  as  follows: 

a.  The  total  volume  of  the  reregulation  pool  is  found  by 

VF  =  VI  +  «ln  -  Qout)At  (28) 

where 

V„  =  final  volume 
F 

V^.  =  initial  volume 

Qin  =  input  to  reregulation  pool  (generation) 

Qout  ~  outflow  from  the  reregulation  pool 

(pumpback  and/or  downstream  release) 

At  =  time  step  length 

b.  From  this  total  volume,  the  elevation  of  the  pool  is 
interpolated  from  a  table  based  on  the  bathymetry  of  the 
pool . 

c.  The  change  in  storage  for  each  longitudinal  segment 
along  the  pool  is  found  through  a  table  relating  the 
pool  elevation  with  volume  for  each  segment. 


d.  The  discharge  is  then  calculated  at  each  interior  cross 
section  bv 


Vi  ■  -  fr 


(29) 


where 

i  =  cross-section  number;  1  =  0  to  number  of 
segments 

=  upstream  discharge  of  segment  i 
AS^  *  change  in  segment  i  storage 

59.  The  advantage  of  this  method  over  the  sump  method  is  that  it 
allows  longitudinal  variation  in  reregulation  pool  temperature  to  be 
modeled.  The  temperature  and  constituent  transport  is  computed  by  using 
this  discharge  distribution  in  an  implicit  finite  difference  transport 
code  to  be  discussed  later. 

60.  The  inclusion  of  level  pool  routing  exclusive  of  other  hydro- 
logic  routing  methods  was  due  to  its  simplicity  and  its  applicability  to 
reregulation  pools.  It  is  noted,  however,  that  the  resolution  of  the 
hydrodynamics  concerned  with  the  propagation  of  the  wave  surge  through 
the  pool  for  operation  changes  cannot  be  determined  with  this  method. 
However,  this  may  not  be  necessary  to  get  reasonable  constituent  and 
temperature  distributions  for  many  problems  (such  as  investigation  of 
hourly  or  daily  average  conditions)  associated  with  reregulation  pool 
hydrodynamics . 

61.  The  level  pool  assumption  may  be  too  severe  if  the  width  of 
the  pool  changes  abruptly  with  water-surface  elevation.  This  would  be  a 
problem  if  the  water-surface  slope  is  great  enough  to  result  in  signifi¬ 
cant  differences  in  surface  area  and  hence  heat  exchange  along  the  longi¬ 
tudinal  axis  of  the  pool.  Impoundments  with  low  velocities  would  have 
very  slight  surface  slopes  making  the  assumption  of  a  level  pool  quite 
reasonable.  As  an  example,  a  uniform  channel  with  a  velocity  of  1  fps 
and  a  hydraulic  radius  of  20  ft  would  have  a  friction  slope  of  about 

7.5  *  10  ^  for  a  relatively  high  Manning's  n  of  0.03.  This  is  less 
than  0.04  ft  per  mile.  Thus,  even  at  this  velocity  and  Manning's  n  , 


the  mildness  of  the  computed  friction  slope  adds  validity  to  the  level 
pool  assumption. 

62.  Another  alternative  involves  solution  of  the  Saint-Venant 
equations  for  the  velocity  currents  and  depths.  For  completeness,  a 
derivation  of  these  equations  similar  to  that  of  Fread  (1983)  is 
included  in  Appendix  A.  In  conservation  form  the  continuity  equation 


|2  +  |A  -  0 

ax  at 


The  momentum  equation  is: 


aq  ,  3(q?/a)  .  ah  .  „ 
at  +  ax  +  gA  a^  +  sf 


where 


2n  n 
c  =  n  Q  Q I 

f  2.21A2R4'3 


A  *  cross-sectional  area 
g  *  acceleration  of  gravity 
h  =  water-surface  elevation 
Sj.  «  friction  slope 
Q  =  discharge 
R  ■  hydraulic  radius 
x  -  distance  along  the  waterway 


n  *  Manning's  n 

63.  The  dynamic  routing  method  is  the  most  rigorous  of  the 
methods  presented.  Water-surface  elevations  and  discharge  at  each  cross 
section  are  calculated  from  the  continuity  and  momentum  equations  of  the 
Saint-Venant  equations. 

64.  The  numerical  hydrodynamic  code  developed  for  this  investiga¬ 
tion  is  based  upon  the  modeling  work  of  Fread  (1978,  1980),  The  model 


is  an  implicit  finite  difference  solution  of  the  Saint-Venant  equations. 
The  development  of  the  finite  difference  formulation  uses  a  four-point 
method  in  which  a  time-weighted  spatial  derivative  and  a  space-centered 
temporal  derivative  is  employed.  The  details  of  the  solution  method  are 
presented  in  Appendix  B. 

65.  The  implicitness  of  the  solution  technique  allows  for  the 
relatively  large  time  steps  which  may  be  needed  for  simulating  the  long 
time  periods  associated  with  the  stratification  cycle  of  a  reservoir. 

The  model  accepts  a  discharge  boundary  upstream  (release  from  the  reser¬ 
voir  or  pumpback)  and  a  discharge  or  a  discharge-head  relationship  for 
the  downstream  boundary. 


Alternative  method  for  simulat- 
ing  reregulation  pool  transport 

66.  The  representation  of  the  reregulation  pool  as  a  single 
element  is  a  poor  representation  for  long  narrow  pools.  The  release 
temperatures  are  simulated  as  a  pool  average  which  makes  the  response  to 
climatic  changes  very  slow.  As  an  example,  consider  a  long  narrow 
channel  with  a  rectangular  cross  section.  The  linear  advection  diffu¬ 
sion  equation  with  surface-heat  exchange  can  be  derived  in  the  same 
manner  as  Equation  27  and  written  as 
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(32) 


where 

U  =  constant,  channel  velocity 
Ej  =  longitudinal  dispersion  coefficient 
Tp  =  the  equilibrium  temperature 
D  =  the  channel  depth  of  flow 

A  large  flow  rate  in  a  long  narrow  pool  renders  the  dispersion  term 
insignificant  compared  to  the  advective  term.  For  the  purpose  of 
example,  consider  the  initial  pool  temperature  to  be  the  equilibrium 
temperature  and  the  temperature  of  the  inflow  to  the  pool  to  be  a 


constant  1°  F  above  the  equilibrium.  may  then  be  subtracted  from  T 

and  the  partial  differential  equation  redefined  as: 


8T  8T  KT 

8t  3x  Dye 


where  the  second  order  dispersion  has  been  neglected.  The  solution  to 
this  equation  is 


_K _  x 

Yc  D  U 


T(t)  =  T  (x  -  Ut)e 
o 


where  L  is  the  length  of  the  pool.  To  calculate  the  release 


temperature,  set 


L  BDL  _  V 

“  U  BDU  Q 


where 

B  =  width 
V  =  pool  volume 
Q  =  discharge 

since  the  temperature  at  the  outflow  is  reached  at  this  time  t  after 
initiation.  Tq  at  the  upstream  end  of  the  pool  is  a  constant  1°  F, 
therefore 


T  =  e 


_K _  V 

Y^D  Q 


Choosing, 


=  140  BTU/ft  -D-°F 
=  62.4  lb/ft3 
=  1.0  BTU/lb-°F 


IB 


V  =  10  KAF 
Q  =  10  KAF/D 

results  in  the  solution  to  Equation  32 

T  =  0.899°  F 

This  temperature  at  the  outflow  requires  a  time  of  t  =  V/Q  =  1  day  for 
the  analytic  solution.  The  simulation  as  a  pool  with  no  longitudinal 
resolution  in  a  manner  similar  to  the  sump  method  previously  used  with 
WESTEX  results  in  the  following  solution,  from  Equation  20: 


T(t)  = 


Q  +  A 
yc 

P  / 


Q  ,  ka 

V  yc  V 
P  - 


+  1  e 


11  +  -ZL-)  t 
ri  V/ 


T(t)  =  0.899 [l  -  e~(KI12)t] 


which  results  in  the  following  tabulation: 

t,  days  T 

1  0.603 

2  0.802 

k  0.889 

8  0.899 

*  0.899 

The  release  temperature  predicted  by  the  sump  method  will  therefore 
change  too  slowly  for  fluctuations  in  flow  temperature  and  equilibrium 
temperature.  It  is  apparent  that  an  improvement  of  the  modeling  method 
to  allow  for  longitudinal  resolution  in  the  reregulation  pool  is  needed. 
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67.  The  formulation  of  the  constituent  transport  model  used  in 
conjunction  with  hydrodynamic  information  and  contributed  from  the  level 
pool  routing  or  the  Saint-Venant  formulation  is  based  upon  the 
advection-dif fusion  equation  discussed  in  paragraph  66,  but  written  for 
a  nonrectangular  channel,  i.e., 


3C  ,  1  3(QC)  1  3(AEL  3x)  . 

.  '  t  <>.  “*  *  — 


A  9x 


where 

C  =  constituent  concentration 
t  »  time 

A  =  cross-sectional  area  at  a  given  longitudinal  location 
Q  =  discharge  through  A 
x  =  longitudinal  distance 
Ej  =  longitudinal  dispersion  coefficient 
S  =  sources  and  sinks  for  the  given  constituent 
The  solution  is  provided  by  an  implicit  finite  difference  method  which 
is  first-order  accurate  in  both  time  and  space.  A  donor  cell  arrange¬ 
ment  is  used  for  the  advective  (9QC/3x)  term.  Again,  as  in  the  hydro- 
dynamic  portions  of  the  code,  the  implicit  formulation  provides  less 
restrictive  time  step  criteria  than  an  explicit  methodology. 

68.  The  surface-heat  transfer  process  (the  predominate  source  and 
sink  terms  for  temperature)  is  solved  by  the  same  basic  method  in  the 
reregulation  pool  and  the  reservoir  models  except  for  one  difference. 

In  the  upper  layers,  the  vertical  distribution  of  heat  is  not  simulated 
in  the  reregulation  pool  model.  The  process  uses  the  approach  developed 
by  Edinger,  Brady,  and  Geyer  (1974),  as  discussed  in  paragraphs  27 
through  31  which  deal  with  analytic  analysis,  with  the  exception  that 
the  surface  temperature  is  assigned  to  be  the  temperature  of  the  partic¬ 
ular  longitudinal  segment.  The  computation  of  equilibrium  temperature 
is  based  solely  upon  meteorological  data  as  outlined  by  Edinger, 
Duttweiler,  and  Geyer  (1968).  The  thrust  of  this  effort  is  in  selection 


of  an  appropriate  method  to  supply  the  hydrodynamics  (water  depths  and 
current  velocities),  therefore  less  detail  is  given  here  for  the  trans 
port  scheme.  However,  a  more  complete  documentation  of  this  transport 
formulation  is  given  in  Appendix  C. 
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PART  IV:  NUMERICAL  COMPARISON  OF  METHODS 

Description  of  Tests 

69.  To  determine  the  adequacy  of  the  modeling  methods,  an  example 
of  a  previous  study  at  WES  was  rerun  for  limited  conditions.  The  advan¬ 
tage  of  using  a  previous  study  was  that  calibration  parameters  were 
already  known.  The  study  chosen  was  the  investigation  by  Holland  and 
Dortch  (1982)  of  the  proposed  reregulation  pool  for  Norfork  Lake, 
Arkansas.  The  proposed  reregulation  facility  is  one-dimensional  in 
appearance.  The  afterbay,  as  previously  modeled,  had  a  volume  of 
12,600  acre-ft  and  covered  an  area  of  470  acres  at  maximum  pool.  Nor¬ 
fork  Lake  itself  impounds  1,295,700  acre-ft  and  covers  22,620  acres  at 
the  top  of  the  April-to-September  conservation  pool,  el  554  ft.  National 
Geodetic  Vertical  Datum  (NGVD) .  The  proposed  reregulation  pool  was  to 
be  about  3.5  miles  long  which  resulted  in  the  use  of  a  worst-case  condi¬ 
tion  of  a  plug  flow  formulation  in  the  sump  method  in  this  and  the  pre¬ 
vious  study.  The  bathymetry  used  in  this  investigation  is  shown  in 
Figure  8. 

70.  The  proposed  project  was  to  be  used  for  pumped-storage ,  so 
the  impacts  of  the  various  methods  upon  the  simulation  of  stratification 
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Figure  8.  Assumed  bathymetry  for  Norfork  Reregulation  Pool 


46 


of  the  upstream  reservoir  as  we]]  as  the  release  from  the  afterhay  may 
be  compared.  Figure  9  reviews  the  schedule  of  inflow,  release,  genera¬ 
tion,  and  pumpback  rates  that  were  simulated.  The  pumpbark  occurred 
during  the  stratification  season,  beginning  about  day  150  (30  April).  A 
description  of  the  modeling  procedure  for  entrainment  and  parameters  for 
modeling  the  upper  reservoir  are  given  in  Holland  and  Dortch  (1982). 

71.  The  one-dimensional  methods  utilized  a  discretization  of 
10  segments.  The  Manning's  n  roughness  criteria,  used  in  the  Saint- 
Venant  equations,  was  set  at  0.03.  The  time  step  used  for  both  the 
Saint-Venant  and  level  pool  methods  was  chosen  as  1  hr.  A  separate  run 
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Figure  P.  Hvdraulir  flow  boundary  information  for  1975  representa¬ 
tive  condition  for  Norfork  Lake 
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using  the  Saint-Venant  method  was  made  using  a  time  step  of  15  min.  The 
results  were  nearly  identical  to  those  using  the  1-hr  time  step  length. 

So,  all  subsequent  comparisons  were  made  with  1-hr  time  steps.  The 

2 

dispersion  coefficient,  ,  was  estimated  to  be  about  2  ft  /sec  based 

upon  the  work  of  Elder  (1959).  The  remaining  parameters  for  entrain¬ 
ment,  heat  exchange,  etc.  were  the  same  for  all  three  schemes  and  were 
identical  to  those  developed  and  used  by  Holland  and  Dortch  (1982)  dur¬ 
ing  the  previous  study. 


Results 


72.  The  results  of  the  simulation  of  reregulation  pool  release 
temperatures  for  each  of  the  three  methods  are  shown  in  Figure  10.  The 
sump  method  showed  average  daily  release  temperatures  to  vary  from  about 
7°  C  in  late  winter  to  almost  18.5°  C  in  the  fall.  The  equilibrium 
temperatures  are  periodically  registered  on  the  plot  as  well.  The  sump 
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Figure  10.  Release  temperatures  from  the  proposed  Norfork 

Reregulation  Pool 


method  mixes  the  net  volume  of  flow  generated  upstream  with  the  reregu¬ 
lation  pool  so  that  the  entire  pool  is  at  one  temperature.  With  this 
method,  the  sump  should  respond  very  slowly  to  temperature  fluctuations. 
For  the  other  two  methods,  the  reregulation  pool  can  develop  longitu¬ 
dinal  gradients  so  that  the  water  leaving  the  system  during  warm  weather 
will  have  had  more  time  to  heat  (or  cool  depending  on  the  time  of  year) 
than  would  water  that  had  just  reached  the  center  of  the  pool.  There¬ 
fore,  the  longitudinal  segmenting  allows  the  model  release  temperatures 
to  be  more  strongly  influenced  hy  the  environmental  conditions.  How¬ 
ever,  the  sump  method  in  this  figure  in  fact  responds  quite  quickly  to 
changes  in  equilibrium  temperature.  It  compares  well  with  the  other 
methods  for  this  case.  The  reason  that  the  comparison  is  good  for  the 
sump  method  here  is  two-fold:  first,  the  reregulation  pool  is  rela¬ 
tively  small  and  the  detention  time  is  only  about  1  day,  and  second,  the 
use  of  the  pumphack  coefficient  has  a  strong  influence  on  the  release 
temperature.  Low  pumphack  coefficients  mix  the  bulk  of  the  generation 
volume  with  the  reregulation  pool  so  the  release  temperatures  will 
respond  slowly  to  climate  changes.  Very  high  pumpback  coefficients 
return  the  generation  flow,  thus  providing  less  mixing  with  the  reregu¬ 
lation  pool  and  allowing  the  predicted  release  temperatures  to  change 
quickly  with  climate  change.  Lnf or  tuna  to  ly  ,  the  pumphack  coefficients, 
while  improving  the  temperature  estimate  for  release,  could  in  fact  pro¬ 
vide  poor  temperature  prediction  in  the  reregulation  pool. 

73.  Two  of  the  main  reservoir  temperature  profiles  computed  dur¬ 
ing  the  period  of  the  greatest  stratification,  July  9  and  September  18, 
are  shown  in  Figure  11.  The  predicted  temperature  profile  for  July  9 
using  the  sump  method  for  the  afterbav  produced  temperatures  between  9° 
and  29°  C.  The  level  pool  results  produced  a  temperature  profile  quite 
close  to  that  computed  with  the  sump  method.  For  this  reason,  the  level 
pool  and  the  Saint-Venant  profile  results  were  plotted  as  deviations 
from  the  sump  method  profile.  The  middepth  layers  of  the  upper  reser¬ 
voir  for  both  one-dimensional  methods  for  modeling  the  reregulation  pool 
were  not  as  warm  as  predicted  by  the  sump  method.  The  difference  was  as 
much  as  0.3°  C.  The  Saint-Venant  method  resulted  in  a  slightly  colder 
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Figure  11.  Upper  reservoir  temperature,  Norfork 
lake  results 

profile  than  the  level  pool  method  by  roughly  0.05°  C. 

74.  The  September  18  profiles  continued  in  the  same  manner  with 
the  one-dimensional  methods  predicting  more  stratification  in  the  main 
reservoir  than  the  sump  method.  The  sump  method  produced  a  temperature 


i  r'  2 


profile  in  the  upper  reservoir  which  ranged  between  13°  and  22.5°  C. 

The  one-dimensional  methods  produced  predicted  hypolimnion  temperatures 
which  were  as  much  as  0.5°  C  cooler.  The  two  one-dimensional  methods 
(level  pool  and  Saint-Venant)  predicted  nearly  identical  profiles. 

75.  It  seems  apparent  that  the  sump  method,  which  has  only  one 
compartment  to  represent  the  reregulation  pool  and  can  only  distinguish 
that  day's  generation  water  quality,  is  lacking  in  many  instances. 
Further,  the  pumpback  coefficient  which  is  critical  to  this  method  is 
somewhat  tenuous  to  determine.  The  pumpback  temperatures  predicted  on 
July  9  were  11.4°  and  10.7°  C  for  the  sump  and  the  Saint-Venant  simula¬ 
tions,  respectively.  The  pumpback  temperatures  predicted  for  Septem¬ 
ber  18  were  15.4°  and  15.0°  C  for  these  two  methods.  The  difference 
could  be  reduced  by  increasing  the  pumpback  coefficient.  This  would 
also  make  the  comparison  between  reregulation  pool  release  temperatures 
and  main  reservoir  strati f ication  patterns  closer.  The  one-dimensional 
methods  have  the  added  flexibility  of  many  segments  to  represent  the 
reregulation  pool.  The  primary  difference  between  the  one-dimensional 
methods  results  from  the  friction  slope  causing  different  water-level 
predictions  in  the  two  methods.  This  water-level  difference  can  allow 
different  surface-heat  exchange  due  to  a  difference  in  surface  area. 

With  the  relatively  short  reregulation  pool  in  this  example  (3.5  miles, 
the  majority  of  the  reregulation  pools  in  Table  1  are  shorter  than  the 
distance),  the  water-level  differences  were  not  substantial  and  the 
results  were  nearly  identical  for  the  level  pool  and  the  Saint-Venant 
methods . 

76.  Tn  order  to  investigate  potential  differences  in  the  one¬ 
dimensional  methods,  the  length  of  the  reregulation  pool  was  increased 
by  five  times  (17.5  miles)  and  the  1-vear  simulation  rerun.  This  change 
should  indicate  the  magnitude  of  the  difference  to  be  expected  between 
results  from  the  level  pool  and  the  Saint-Venant  method  results.  The 
upper  reservoir  temperature  profi’es  for  lulv  9  and  September  18  dates 
are  shown  in  Figure  12.  The  temperatures  of  the  metal imnion  and  hvpo- 
limnion  for  both  dates  were  as  much  as  0.6°  f  warmer  as  predicted  with 
hydrodynamics  supplied  bv  the  Saint-Venant  method  compared  with  level 
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Figure  12,  Temperature  profiles  in  the  main  reser¬ 
voir  for  a  17.5-mile-long  reregulation  pool 

pool  routing.  The  release  temperatures  for  these  two  methods  are  shown 
in  Figure  13.  Again,  the  temperatures  predicted  bv  the  Saint-Venant 
method  were  warmer  bv  about  1°  C.  The  overall  comparison  between  the 
two  methods  was  quite  close  with  over  99  percent  of  the  variance 


explained.  The  differences  in  the  relative  temperatures  from  the  level 
pool  to  Saint-Venant  simulations  are  a  function  of  channel  slope,  the 
channel  cross-section  geometry,  and  the  channel  roughness.  For  the  same 
volume,  the  average  depth  may  be  increased  or  decreased  by  the  level 
pool  assumption  relative  to  the  Saint-Venant  results.  This,  in  turn, 
affects  the  rate  of  heat  exchange  through  the  water  surface.  Further, 
the  Manning's  n  roughness  factor  used  in  the  Saint-Venant  simulation 
is  a  calibration  parameter  which  is  estimated.  In  this  case,  the  simu¬ 
lation  used  a  value  of  0.030  which  is  probably  somewhat  high  for  a 
straight  prismatic  channel.  As  a  sensitivity  check,  the  year  simulation 
was  rerun  with  a  Manning's  n  value  of  0.015.  The  results  are  included 
in  Figure  13  for  release  temperatures.  These  results  were  quite  close 
< the  level  pool  method  was  less  than  0.2°  C  lower  on  average)  to  those 
of  level  pool  routing,  even  for  this  relatively  long,  shallow  and  narrow 


reregulation  pool, 


.sW  \  h 

Ct  \  i  \ 


LEGEND 


-  IE VE l  POOL 


•  •  •  ST  VENAMT  In  0  030) 

&  &  &  ST  VENANT  In  00151 

O  EQUILIBRIUM  TEMPERATURE 


?00 

CALENDAR  DAY 


Figure  13.  Release  temperature  from  the 
17.5-mile-long  reregulation  pool 
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Discussion 


77.  The  Saint-Venant  simulations  often  proved  to  be  difficult  to 
run.  These  difficulties  arose  as  a  direct  result  of  boundary  conditions 
which  are  typical  of  reservoir/rerc gulation  pool  modeling.  Instead  of 
the  usual  riverine-type  boundary  conditions,  in  which  the  discharge  is 
given  upstream  and  water  levels  or  a  stage-discharge  relationship  down¬ 
stream,  discharge  is  most  often  known  both  upstream  and  downstream  with 
no  water-level  information  known  at  either  location  for  the  case  of 
interest.  This  causes  problems  in  conserving  mass,  due  to  computational 
roundoff  and  precision.  This  problem  is  magnified  by  the  long  simula¬ 
tion  periods  required.  The  problem  may  be  demonstrated  by  considering  a 
situation  in  which  a  discharge  Q  enters  the  basin  for  12  hr  (with  no 
reregulation  outflow)  and  leaves  the  basin  (with  no  additional  inflow) 
at  the  same  rate  for  the  next  12  hr.  The  water  volume  at  the  beginning 
and  end  of  the  day  should  be  identical.  However,  as  one  computationally 
steps  through  the  day,  calculating  the  discharge  distribution  in  the 
basin  and  the  water  elevation  by  the  equations  of  continuity  and  momen¬ 
tum  and  the  bathymetric  tables,  the  volume  at  the  end  of  the  day  will 
not  precisely  match  that  at  the  day's  beginning.  After  a  year  of  simu¬ 
lation,  this  could  amount  to  20  ft  or  more.  The  long  periods  used  in 
reservoir  investigations  cause  this  to  be  a  significant  problem.  This 
problem  was  alleviated  to  some  degree  by  a  small  convergence  criterion 
(flow  rate  within  0.1  percent,  and  water-surface  elevation  less  than 
0.005  ft)  and  additional  coding  which  checks  the  basin  volume  after  each 
operation  cycle  and  makes  small  flow  corrections  to  prod  the  calcula¬ 
tions  toward  more  precise  conservation  of  mass. 

78.  The  level  pool  method  gave  results  which  were  comparable  to 
the  more  elaborate  Saint-Venant  results,  even  for  the  relatively  long 
narrow  case  of  this  reregulation  pool.  The  WESTEX  code  using  the  level 
pool  routing  required  only  a  third  the  processor  time  that  the  WESTEX 
code  using  the  Saint-Venant  scenario  required.  The  Saint-Venant  solu¬ 
tion  does  not  lend  itself  well  to  the  particular  boundary  conditions 
associated  with  reservoir  simulations. 


79.  The  three  methods  considered  for  simulation  of  hydrodynamics 
which  drive  a  transport  simulation  were  the  sump  method,  level  pool 
routing,  and  modeling  using  the  finite  difference  representation  of  the 
Saint-Venant  equations.  These  methods  were  compared  by  simulation  of  a 
previous  study  of  the  proposed  reregulation  pool  at  Norfork  Reservoir, 
Arkansas.  As  the  proposed  facility  was  to  provide  pumped-storage  capa¬ 
bility  for  the  project,  release  temperatures  as  well  as  the  vertical 
temperature  distribution  in  the  main  reservoir  were  compared  for  each 
modeling  method.  The  reregulation  pool  simulation  methods  were  coupled 
with  the  vertical  one-dimensional  thermal  model  WESTEX  to  simulate  the 
upper  reservoir. 

80.  An  analytical  comparison  of  the  sump  method  with  the  other 
two  methods  demonstrated  that  the  lack  of  resolution  in  the  reregulation 
pool  with  the  sump  method  resulted  in  a  slow  response  time  for  predicted 
release  temperatures.  If  pumpback  is  being  used  at  the  site,  a  pumpback 
coefficient  improves  the  prediction  of  release  temperatures.  The  effect 
of  pumpback,  as  predicted  by  the  sump  method,  was  less  stratification  in 
the  upper  reservoir  during  the  stratification  season  than  that  predicted 
by  the  other  two  methods.  The  ability  of  the  sump  method  to  produce 
reliable  results  depends  upon  the  selection  of  a  pumpback  coefficient 
that  generally  requires  field  or  physical  model  input.  The  need  for 
this  coefficient  is  again  caused  by  the  lack  of  resolution  of  the  sump 
method.  This  pumpback  coefficient  determines  the  amount  of  the  previous 
generation  flow  for  a  particular  day  which  is  returned  to  the  main  reser¬ 
voir  during  the  day's  pumpback.  With  proper  resolution  this  would  be 
determined  automatically. 

81.  The  length  of  the  pool  for  which  these  comparisons  were  made 
was  3.5  miles.  However,  even  for  this  relatively  long,  narrow  pool, 
there  was  essentially  no  difference  between  the  predictions  for  the 
vertical  upper  reservoir  temperature  profiles  and  the  release  tempera¬ 
tures  from  the  reregulation  pool  made  using  level  pool  routing  and  the 
full  Saint-Venant  equations.  A  hypothetical  reregulation  pool 

17.5  miles  long,  with  the  same  widths  and  cross-sectional  areas  as  pre¬ 
viously  modeled,  was  simulated  with  these  two  methods.  The  results 


showed  some  differences  between  predicted  release  temperatures  for  the 
two  one-dimensional  methods.  However,  these  differences  were  so  small 
as  to  be  masked  by  the  estimate  of  Manning's  n  in  the  Saint-Venant 
sirml a t i ons . 

82.  The  finite  difference  formulation  of  the  Saint-Venant  equa¬ 
tions  also  suffered  difficulties  in  the  form  of  lack  of  conservation  of 
mass  with  the  specified  upstream  and  downstream  disci, urge  boundary  con¬ 
ditions  typical  of  reservoir  modeling.  The  conservation  of  mass  prob¬ 
lems  were  due  to  the  imprecision  of  calculaf ions  and  interpolation  of 
the  elevation /volume  tables  and  the  long  time  periods  of  simulation, 
resulting  in  poor  conservation  of  mass.  Additional  coding  reduced  the 
problem  by  making  small  corrections  to  the  given  flow  rates  to  account 
for  the  error  in  volume  in  the  reregulation  pool . 
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PART  V:  SUMMARY  AND  CONCLUSIONS 


Summary 

83.  A  reregulation  dam  is  used  to  moderate  flow  fluctuations  due 
to  reservoir  releases,  often  as  a  result  of  hydropower  operations.  The 
pool  formed  by  the  reregulation  dam  can  also  serve  as  the  water  volume 
used  for  pumpback  in  a  pumped-storage  generation  system.  The  reregula¬ 
tion  pool  will  impact  the  temperature  of  the  reservoir-reregulation  pool 
system  and  the  downstream  releases.  A  cost-effective  method  for  simu¬ 
lating  the  hydrodynamics  and  heat  transport  were  presented.  The  method 
was  coupled  with  a  vertical  one-dimensional  reservoir  model  to  simulate 
the  hydrodynamics  and  heat  transport  of  the  reservoir-reregulation  pool 
system. 

84.  The  modeling  procedure  for  simulating  the  hydrodynamics  and 
heat  transport  of  the  reregulation  pool  was  selected  from  three  candi¬ 
date  methods.  The  three  methods  in  progressive  order  of  complexity  are: 
(a)  the  sump  method,  (b)  the  level  pool  routing  method,  and  (c)  the 
finite  difference  simulation  of  the  Saint-Venant  equations.  The  methods 
were  incorporated  within  the  WESTEX  model  and  the  predicted  temperature 
profiles  in  the  main  reservoir,  and  the  release  temperatures  from  the 
reregulation  pool  were  compared  to  results  from  the  Norfork  Lake  study 
previously  conducted  at  WES. 


Conclusions 


85.  As  demonstrated  by  the  analytic  descriptions,  the  resulting 
conclusions  are  given  below. 

a.  The  increased  detention  time  provided  by  a  reregulation 
pool  will  result  in  an  increased  warm-season  temperature 
in  the  reach  In  which  the  pool  is  placed.  The  pool  is 
generally  sufficiently  deep  that  short-period  fluctua¬ 
tions,  such  as  the  diurnal  temperature  range,  will  not  be 
significant . 


b.  The  more  uniform  flow  leaving  the  reregulation  pool  will 
result  in  a  moderation  of  water  temperature  fluctuations 
downstream  of  the  reregulation  dam,  compared  with  condi¬ 
tions  without  the  pool. 

£.  A  pumped-storage  system  causes  a  two-way  interaction 
between  the  major  reservoir  and  the  reregulation  pool. 

An  analysis  of  a  homogeneous  system  demonstrates  that  the 
similarity  of  conditions  in  the  upper  reservoir  and  the 
reregulation  pool  increases  with  increasing  amounts  of 
pumpback.  A  conservative  constituent  found  in  the  main 
reservoir  will  reach  its  peak  concentration  in  the 
reregulation  pool  sooner  and  with  a  higher  concentration 
for  increasing  generation/pumpback  flows  between  two 
homogeneous  reservoirs. 

d.  The  entrainment  and  mixing  processes  associated  with 
pumpback  will  result  in  a  more  uniform  density  profile  in 
the  upper  reservoir  than  that  observed  without  pumpback. 
The  degree  of  influence  of  the  entrainment  and  mixing  is 
dependent  upon  the  relative  size  of  the  reservoir,  the 
rate  of  pumpback  flow,  and  the  density  structure  of  the 
upper  reservoir. 

e.  The  complexity  of  geometric  description  and  the  physical 
processes  makes  modeling  necessary  when  investigating  the 
effects  of  a  reservoir-reregulation  pool  system.  The 
time  period  necessary  for  this  simulation  is  that  of  the 
natural  stratification  cycle,  which  is  generally  a  year. 
An  efficient  simulation  routine  for  such  a  long  period 
should  be  used  to  minimize  computational  costs. 

£.  The  very  large  detention  times  in  the  upper  reservoir  and 
the  small  detention  times  of  the  existing  reregulation 
pools,  as  well  as  historical  modeling  practice,  suggest 
that  resolution  of  vertical  stratification  patterns  is 
critical  in  the  upper  reservoir  while  longitudinal  gradi¬ 
ents  are  of  most  importance  in  the  reregulation  pool. 

86.  The  results  of  these  simulations  indicated  that  the  level 
pool  method  was  generally  the  most  efficient  modeling  method  of  the 
three.  It  provided  sufficient  resolution,  was  easier  to  run,  was  less 
expensive,  and  supplied  results  comparable  to  the  more  sophisticated 
Saint-Venant  method.  The  sump  method  requires  the  determination  of  a 
pumpback  coefficient  to  provide  good  comparison  of  release  temperatures. 
In  the  event  that  pumpback  is  not  used,  the  prediction  of  release 
temperatures  will  change  slowly.  The  Saint-Venant  method  appears  to  be 
needed  only  if  there  is  a  substantial  water-surface  slope,  if  the  flow 
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Carters  Reregulation  Pool 
Maximum  storage  17.5  KAF 
Surface  area  at  maximum 
storage  of  0.87  KA 
Length  1.5  miles 

DeGray  Reregulation  Pool 
Maximum  storage  3  KAF 
Surface  area  at  maximum 
storage  0.43  KA 
Length  2  miles 

Clarence  Cannon  Reregulation  Pool 
Maximum  storage  8.4  KAF 
Surface  area  at  maximum 
storage  1.02  KA 
Length  9.5  miles 

B.  A.  Steinhagen  Lake 
Storage  at  interim  level  94.2  KAF 
Surface  area  at  interim 
level  13.7  KA 
Length  30  miles 

Broken  Bow  Reregulation  Pool 
Design  pool  storage  9.5  KAF 
Design  pool  surface  area  0.53  KA 
Length  3  miles 

Keystone  Reregulation  Pool 
Maximum  storage  3.5  KAF 
Length  8  miles 


Bonneville  Lake 

Normal  maximum  storage  537  KAF 
at  normal  maximum  surface 
area  20.4  KA 
Length  4.5  miles 

Foster  Reservoir 
Maximum  storage  61  KAF 
Surface  area  at  maximum 
storage  1.22  KA 

Big  Cliff  Reregulation  Pool 
Total  storage  6.45  KAF 
Surface  area  at  maximum 
pool  0.141  KA 
Length  2.8  miles 

Dexter  Reregulation  Pool 
Total  storage  27.5  KAF 
Surface  area  at  normal 
pool  1.03  KA 
Length  2.7  miles 


Carters  Reservoir 
Maximum  storage  473  KAF 
Surface  area  at  maximum 
storage  3.22  KA 


DeCrav  Lake 

Maximum  storage  655  KAF 
Surface  area  at  maximum 
storage  17  KA 


Mark  Twain  Lake 
Maximum  storage  2,060  KAF 
Surface  area  at  maximum 
storage  38.4  KA 


Sam  Rayburn  Dam  and  Reservoir 
Maximum  storage  4,000  KAF 
Surface  area  at  maximum 
storage  143  KA 


Broken  Bow  Reservoir 
Maximum  storage  1,370  KAF 
Surface  area  at  maximum 
storage  18  KA 

Keystone  Lake 
Maximum  storage  1,840  KAF 
Surface  area  at  maximum 
storage  55.3  KA 


The  Dalles 

Maximum  storage  at  200,000  cfs 
859  KAF 


Green  Peter  Reservoir 
Maximum  storage  430  KAF 
Surface  area  at  maximum 
storage  3.72  KA 

Detroit  Reservoir 
Maximum  storage  455  KAF 


Lookout  Point  Reservoir 
Normal  storage  42C  KAF 
Normal  surface  area  4.2  KA 


Coosawatter  River,  Georgia 
Average  flow  cfa 

Hvdropower  pumped-atorage  svstem 


Caddo  River,  Arkansas 

Average  flow  726  cfs 

Hydropower  pumped-etorage  system 


Salt  River,  Missouri 
Average  flow  1,530  cfs 
Hydropower  pumped-storage  svstem 


Angelina  River,  Texas 
Average  flow  2,800  cfa 
Maintains  flow  to  reduce  sallnltv 
Intrusion 


Mountain  Fork  River,  Oklahoma 
Average  flow  1,300  cfs 
Hvdropower  reregulatior  and  to 
satlsfv  lev  flow  requirements 

Arkansas  River,  Oklahoma 
Average  flow  fc.ifp  <-fs 
Built  to  maintain  minimum  f  1  nv* 
water  quail  tv  at  Tulsa  and  to 
moderate  hvdr  -power  releases 

Columbia  River,  Washington  a'*** 
Oregon 

Average  flow  Q  4,3^0  r  *  « 
Reregulatlor  to  meet  r a  1 
rest  r lc  t 1  or  « 

South  certla-  Q  t  •  :»t  .  rec=  •“ 
Average  • ’  w  ’ , • v  r  *  « 

Provides  f  1  po  ;  n  t  r ’  a »*  1 
power  re  re?  -  '  1 1  1 

North  Santlar  Mver.  rtrr 
Average  flow  , ^ cfs 
uvdropnver  reregul at  * -r 


Middle  F n r k  V 
nreRi — 

Average  f  \  ow  .  V  *■ 
hvdropower  reregu'a 
power  production 


amef  r  •  u  1  v#  r 


NOTE:  The  Dardanelle  Reservoir,  Arkansas,  was  also  mentioned  as  bavin?  an  -  W1 

ever  it  is  rarely  used  in  that  capacity.  KA  •  Klloacre,  KAF  •  K  i a  re  -  ‘  eer  .  '«  • 
par  second. 
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Reregulation  Pool  Site 


Temperature  Amplitude 
(Water  Temperature  Ampl i tude/Equi 1 ibr ium 

_ Temperature  Amplitude) _ 

Cycle  Period,  days  Reregulation  Pool  Natural  Conditions 


APPF.NO  IX  A:  SAINT-VENANT  EQUATIONS 


I.  The  Sa  int -Ven.mt  equations,  in  the  form  used  in  this  study, 
mav  he  derived  as  follows: 


.  r  t  i  on  ,  t  the  stream  with  cell  walls 
riie  1  innel  bottom.  The  downstream  is 
t  ■'  t  i  t  ream  a>-  ' .  The  water  positions 

Ira  i  ne  ....  water  pcs  i  t  I  on  at  section  fill  at 
;  .■  r  .  '  lir.-wise,  and  !  indicate 


*  f  .  ■  f  I  ' 


:  n.i  .mi;  '  e  ■  s  i  h  1  e  fluid,  described 
.  ■•  >  time  perl  od  ‘  t  ,  won  1  <i 
t  1 1  h  r  r  a  •  t  I  n  the  f  1  ow 


■.•me-  ’  r.  nich  riie  downstream 


c  • .  >  1  a  y  I  v  e  n 


wv' 


$■ 


parameter  over  length  Ax  (i.e.,  A  represents  the  average  instanta¬ 
neous  cross-sectional  area  over  length  Ax)  ,  then 

A(t ) Ax  -  A(t  -  At ) Ax  -  (Q  -  Qn)At  (Al) 

u  D 

The  sign  convention  is  chosen  for  x  positive  to  the  right.  Dividing 
by  AxAt  yields 


A  ( t )  -  A  ( t  -  At) 
At 


(  A2 ) 


Noting  that  Q  -  Q  is  the  same  as  Q(x  +  Ax)  -  o(x)  and  assuming  A 
and  Q  are  sufficiently  smooth  over  the  region  Ax  ,  and  Ax  and  At 
approach  0  .  Then 


ii  ♦  4U  •  o 

At  Ax 


(AD 


This  is  the  continuity  equation. 

3.  Conservation  of  momentum  means  that  the  net  rate  ot  momentum 
entering  an  element  plus  the  sum  of  the  forces  acting  on  the  element 
equals  the  rate  of  accumulation  of  momentum. 

4.  Kepresenting  the  net  rate  of  momentum  entering  an  element  on 
the  upstream  face  is 


if  a  uniform  velocity  distribution  is  assumed  then, 


s 


the  momentum  f  lux  through  a  section  is 


Q(pu)  -  q(pS).  p(^) 


5.  Representing  the  new  rate  of-  momentum  entering  an  element  on 
the  downstream  face  is 


-P 


6.  The  forces  acting  on  an  element  are  listed  below, 
a.  Pressure  forces: 


P  =  p  _  p  +  p 

UDT  U  D  rT 

where 

P^t  *  pressure  force  on  the  upstream  face 

Pjj  *  pressure  force  on  the  downstream  face 

P^.  ■  x  component  pressure  forces  on  the 
sidewalls  of  the  element 

b.  Weight  force: 

W  *  gS  AAx  (A- 

o 


This  is  a  component  of  weight 
axis.  So  is  the  bottom  slope 
between  the  x  direction  and 


force  directed  along  the  x 
equal  to  sin  of  the  angle 
the  horizontal. 


Friction  force: 


D  «  -Pg  S^AAx 


(A 


7.  The  momentum  within  an  element  can  be  expressed  as 


M ( t )  -  pQ(t)Ax 


(Ai 


8.  Therefore,  the  discrete  conservation  statement  would  be 


[c\\  ■ 


-  —  1  +  ogS  A  Ax  -  ogSfAAx 
to  t 


+  -  PD  +  PT  I  -  >  Q(  t )  Ax  -  o  Q  ( t  -  At)  Ax  ( A7 


If  the  flow  depths  y  and  widths  B  change  gradually  over  the  length 
of  the  cell,  then  the  derivatives  SB/ r-x  and  <Y /  ix  exist.  Further, 
if  the  pressure  distribution  is  assumed  hydrostatic  and  that  the  bottom 
slope  is  small,  then  the  pressure  force  on  a  face  is  given  bv : 


P  ”  ug  /  (Y  -  n)  B(x,’l  dr 


where  n  is  0  at  the  channel  bottom  and  Y  at  the  surface 


PT“  og.'ix  /(Y  -  "1 
o 


P  -  P  +  p  -  P  -  —  ‘,x 
U  I)  f  T  1  a  x 


R‘‘X  ^  S(* 


from  Liebnitz’s  Rule 


OgAx  --  /  f Y  -  ' ) Bd' 
o 


g  xl/ 


I Y  -  ■  )(B)d- 


♦  i  Y  Y  1  H  i  \  ■  - 


•  • 


5'V' «.•  -/  .•  -V  •  V  •  •  .  •  • 


ox  *  pgAx  f  Bdn  +  f  (Y  -  n)  dn 
x  o  o 


PT  -  ff  Ax  =  ogAx  y  (Y  -  n)  If  dn  -  /  £  Bdr 
o  o 


t  <?  -  -»  H  ■" 

o 


P.J.  -  —  Ax  -  AgAx  /  (Y  -  n)  —  dr'  -  /  ^  Bdr 


-  /  (Y  -  r;>  —  dr 
-*x 


fY  -  (>„  4  >’T  -  *•-, 


*P  ,  f  <*Y  '  .  ~  *Y 

Tx  Ax  -  -CR..X  f  Jx  B  d-  -  -t  g.  xA  - 

o 


'4.  If  the  weight  force  is  included  with  the  pressure  force 


!>,.  -  l'n  4  lT  +  W  -  -  RAXA 


h  is  tin'  w.i  t  <*  r  -  snr  f  ;i<  e  elevation. 

1 1 .  1 >  i  v  I  <1  i  n  v  t  !>**  d  i  s<  r p  t  p  rmimet)  J  uni  ecjiia  t  i  >  >ri  1'  v  .  y  '  x  *  t  Ip. 
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As  Ax  and  At  +  0  ,  and  substituting  the  weight  and  pressure  force 
derivation  and  rearranging  leaves: 

2 

an  a (n  t a1  /ah  \ 


0 


( 
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APPENDIX  B:  HYDRODYNAMIC  CODE 


1.  The  development  of  the  finite  difference  formulation  uses  a 
four-point  method.  The  time  indicator  is  given  in  the  following 
equations  by  a  superscript  "N"  ,  where  N  represents  the  n  time 
step.  Similarly,  distance  is  given  by  the  subscript  "i"  .  The  spatial 
derivatives  for  a  segment  employed  a  time  weighting  factor  0  which  may 
vary  from  0.5  to  1.  A  value  of  0.5  provides  equal  weighting  of  the 
values  at  the  N  and  N  +  1  time  steps.  Thus,  the  finite  difference 
formulation  for  continuity  used  is: 

-  <n  *  <>  -  «>(<$♦,  -  o 


and  the  momentum  equation  representation  is: 


<Sti_  HJL 


-.■VV\ 


-N 

+  ra 


N  lhN 

i  \N 


i+1  ’  h\  + 


(B2) 


whe  re 


<\+,  *  V 


"i  Qi  Qi 


-2  4/3 

2.2IAJR 


bl 


<.v-  ^  .'w  .a^  •  /a 


» 


These  calculations  are  made  for  each  segment  over  the  length  of  the 
simulated  reach  in  combination  with  the  boundary  conditions  from  which 
the  values  of  Ah^  ,  Ah^+j  ,  AQ^  ,  and  AQ^  +  j  may  found  by 

solving  the  system  of  equations  resulting  from  Equations  B5  and  B6.  T 
superscript  N+l  is  understood  throughout  if  no  superscript  is  given. 
The  various  partial  derivatives  are  evaluated  as  follows: 


3C  Ax  8Ai 

3h  ~  At  3h 


Bi+  ah] 


3(-Qi) 


0.5 

At  3h1+1 


°-5£  K+i  +  mT 


5  (“-^7 


+ 5  h‘)(’ 


h1+]  -  hlAx  Sf 


+  +  AX  (B1 


—  .  n. 


->s 

S  +  -2(?)  +  *A,  Ax  — r1 
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3M 


ah 


i+l 


□  +!!j±ih  \  +  »/B  +!!i±l 

Bi+1  3h, . ,  hi+l  +  2  Bi+1  8h 


i+l 


i+l 


i+l 


x  (h1+1-h1  +  AxSfJ+gAill+AxSfi 


(B12) 


— ^ m  0  5  —  +  0 

3Qi+l  At  + 


as. 


2(  a;i+1  +  *Ai  AX  -gq- 


(B13) 


.  .  .  N+l  lN+1  aN+1  AN+1  ,  .  , 

The  new  estimates  for  h^  ,  ,  or  Qi+1  were  found  by 


f^  (new)  «  f^  (old)  +  Af 


(B 1  4) 


,  ,  ,  ,  ,  N+l  ,  N+l  _N+1 

where  f  is  arv  one  of  h,  ,  h,,,  ,  Q 


,,N+ 1 

or  values. 


If 


i  ■  "i+l  ’ 

f  is  sufficiently  small,  the  f  values  are  specified  as  solutions. 

If  convergence  has  not  been  reached,  these  new  estimates  are  used  to 
generate  updated  estimates  with  the  procedure  just  described. 

4B 

1.  The  terms  involving  —  h  in  this  method  were  original lv 

c*  n 

neglected;  however,  it  was  quickly  discovered  that  changes  in  depth  and 
width  in  a  reregulation  pool  occur  much  too  rapidly  to  drop  these  terms. 

.  To  close  the  svstem  of  equations,  conditions  on  external 
boundaries  must  be  specified.  In  this  case,  a  discharge  historv  Q(tl 
was  designated  both  upstream  and  downstream.  This  condition  can  be 
handled  as  demons t rat ed  bv  Fread  (1483)*  with  the  boundarv  equation 


qJ+1  -  Q(ti 


i  8!  SI 


References  cited  in  the  appendixes  are  included  in  the  reference-,  at 
the  end  ot  the  main  text. 


FU 


where 


9G 

3Q] 

3G_ 

3h, 


1 

0 


for  the  upstream  boundary  (i  *  1)  and 


ac 

3Q 


L 


1 


3C 


0 


for  the  downstream  boundary  (i  «  1). 

5.  The  coefficient  matrix  is  handed  with  most  of  the  elements 
being  zero  except  for  four  elements  in  each  row  ad4acer.t  ti  the  m a;n 
diagonal.  An  efficient  solution  technique  developed  bv  Tread  ' ^ ' 1 > 
which  the  computations  do  not  involve  zero  elements  is  used.  This 
reduces  the  number  of  matrix  operations  '  rom  .  -  *  *»\  *  -  ;N 

i  1 HN  -  ;  0  '■  wiie  re  S  is  the  t  't  a  1  "umber  or  •  r  sect  i  'ns  i '  -ug  r L-  e 

rod:  further  .  the  code  <t  •:  es  -nlv  t  •  e  nor/e*-  e  erne  ts  and  s  »ed\; 


APPENDIX  C:  TRANSPORT  CODE 


' .  The  partial  differential  equation  approximated  by  this  code 
for  temperature  is 


MAT)  j(QT)  3  3 (AT)  “  T) 

;t  3x  ?x  LL  3x  "  Yc 

P 


(Cl) 


where 

A  *  cross-sect ional  area 
T  ■  temperature 
C  *  discharge  rate 

F.  *  longitudinal  dispersion  coefficient 
K  =  surface-heat  exchange  coefficient 
“  =  top  width 

=  equilibrium  temperature 
■  *  weight  densitv 
*  specific  heat 

'.r  r  •  »  •  i  r  i  t  e  difference  formulation,  the  spatial  derivatives  are 
-e:-  t  e:  hv  entered  differences  in  the  dispersion  term  by  Roache's 
:  r  s  t -  •  r  d  t -  windward  differencing  scheme  for  the  advective  term.  The 
w : v dw i r i  scheme  preserves  the  transportive  property  but  is  only  of  first 
tier.  The  temporal  derivative  is  also  only  of  first  order.  So  the 
■•vera  1  scheme  is  approximately  (At, Ax).  The  scheme  used  is  implicit  so 
t^at  the  time  step  restriction  Av/Q  is  removed  where  Av  <=  incremental 
segment  volume.  The  finite  difference  formulation  therefore  is: 


0.5  (A^  + 


where  - 
and  -  0.5  (B±  +  B1_1) 

The  superscript  N  indicates  the  specific  time  step  and  i  indicates 
the  longitudinal  location  as  shown  below. 


2.  All  information  such  as  Q  ,  A  ,  and  V  ,  which  are  supplied 
by  the  hydrodynamic  model,  are  defined  at  the  segment  ends  as  shown  in 
the  figure  for  the  Q  parameter.  The  constituent  concentration  or 
temperature  is  defined  at  the  center  of  the  segment. 

3.  All  information  used  in  this  formulation  is  known  from  the 

hvdrodvnamic  model  results  and  from  previously  calculated  temperatures 
except  for  Tj  .  ,  ,  and  for  each  segment  i  .  The 

'efficients  of  these  parameters  are  then  grouped  so  that  the  resulting 
■« v s t em  of  equations  has  the  form: 


APPENDIX  D:  NOTATION 
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A 

n 

B 

c 

P 

C 

D 

D 

v 

e 


E 

f 

F 


L 

i 

D 

g 

h 


H 

s 

i 

I 
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K  and  T^ 
KA 
KAF 
KAF/D 
L 
n 

0 

P 


Cross-sectional  area;  horizontal  area  of  the  element, 

2 

ft  ;  cross-sectional  area  at  a  given  longitudinal 
location 

Surface  area  of  reservoir  n 
Width 

Specific  heat  of  water 
Constituent  concentration 
Channel  depth  of  flow 
Vertical  entrance  dilution 

Average  normalized  density  gradient  taken  as  0.3 
x  10~6/ft 

Longitudinal  dispersion  coefficient 

Interfacial  friction  factor 

Densimetric-Froude  number 

Gravitational  acceleration  constant 

Water-surface  elevation  above  datum;  water-surface 
elevation 

Average  pond  depth 
Net  rate  of  surface-heat 
Cross-section  number 
Inflow 

Coefficient  of  surface-heat  exchange;  vertical  diffusion 
2 

coefficient,  ft  / day 
Surface-heat  exchange  component 
Kiloacre 
Kiloacre-f eet 
Kiloacre-f eet/D 
Total  pond  length 

1  or  2,  denoting  the  main  reservoir  and  reregulation 
pool,  respectively;  Manning's  n 

Outflow 

Pond  number 


Pressure  force  on  the  downstream  lace 


l' 

Q 


Ni 

Qin 

Qo 

^out 

Qv 

Ql 

Qo 


R  <  R 


c 

S 


X  component  pressure  forces  on  the  sidewaUs  of  the 
element 

Pressure  force  on  the  upstream  face 

Flow-through  rate  In  cfs;  flow  through  the  reach; 
discharge 

Inflow  rate,  cfs/day;  upstream  discharge  of  segment  1 

Input  to  reregulatlon  pool  (generation) 

Condenser  flow  rate;  withdrawal  rate,  cfs/dav 

Outflow  from  the  reregulation  pool  (pumpback  and/or 
downstream  release) 

Net  vertical  flow  into  or  out  of  a  layer 

Inflow  and  outflow  rate;  flow  rate  in  generation  (from 
reservoir  1  to  2) 

Pumpback  rate  (flow  rate  from  2  to  1) 

Net  flow  rate  through  entire  system  (equal  to  the 
inflow  to  reservoir,  and  also  the  release  from 
reservoir  2) 

Ratio  of  yearly  reservoir  inflow  volume  to  actual 

reservoir  volume  for  a  number  of  impoundments  in  the 
continental  United  States;  hydraulic  radius 

Criterion  for  presence  of  reservoir  stratification 

Storage;  sources  and  sinks  for  the  given  constituent 

Friction  slope 

Source  or  sink 

Time 

Water-surface  temperature;  temperature 
Cycle  period 

Intake  temperature  (from  pond);  temperature  of  the 
inflow,  °F 

Temperature  of  reservoir 

Discharge  temperature  (into  pond);  temperature  of  the 
withdrawal,  °F 

Equilibrium  temperature 

Inflow  temperature 

Instantaneous  temperature  of  the  reservoir  and  also  the 
outflow  temperature 


,-vJ 


D2 


Initial  t empe r n r n r e  of  reservoir  I 

I  Itlmate  s  t  (“iii!  v  -  s  t  at  e  t  priper*tun>  for  train  reservoir 
'nltlal  temperat  ore  of  reservoir 

1' Itlmate  st  eadv-staf  e  temperature  of  reregit  1  at  Ion  po 
Constant,  rhanne !  velocity 

Volune  <>f  the  reservoir  In  ruM<  feet;  velocity;  poo 
volume;  volume  of  reach 

Volume  of  reservoir  n 

Finn?  volume 

Inf flal  volume 

Volume  n*  the  reservoir 

Average  pond  width 

Uelghtir.,  roe  f  f  i  r  i  ent  ;  dlatanre  along  the  waterway; 
longitudinal  distant  e 

Kqul  librium  temperature  as  a  ftmrtinn  of  time 
Coefficient  of  thermal  expansion  for  water 
Specific  weight  of  wnfer 

Change  of  storage  within  the  reach  during  a  time 
Inc  rement 

Change  In  segment  1  storage 
Time  step  length 

Total  temperature  difference  IT  -  T  >  across  pond 
Cycle  frequency  -  ? ’  'T 

"i 

Density  of  water,  lh/ft 
Distance  along  longitudinal  axis 


